AFAL  TR-87-043 


AD: 


rTm  f-VT^TT  1T"VT 


m  HLE.coea 


X.li  ■  IPW 

§ 


AD-A18S  808 


Final  Report 
for  the  period 

30  Sept  1985  to 

31  March  1987 


September  1987 


Crack  Opening  and  Extension  in 
Inert  Solid  Propellant  _ 

An  Experimental  Study  _  J- f 

SLECTh 

JAN  0  5  1988 


$ 


Authors: 

D.  Post 
C.  W.  Smith 
in  Collaboration 
with 

R.  Czarnek 


Virginia  Polytechnic  Institute  &  State  Univ. 
Dept,  of  Engineering  Science  &  Mechanics 
Blacksburg,  VA  24061 


F0461 1 -85-C-0098 
352-236-1 


Approved  for  Public  Release 


Distribution  is  unlimited.  The  AFAL  Technical  Services  Office  has  reviewed  this  report,  and  it  is 
releasable  to  the  National  Technical  Information  Service,  where  It  will  be  available  to  the  general 
public,  including  foreign  nationals. 


prepared  for  the:  Air  Force 

Astronautics 

Laboratory 

Air  Force  Space  Technology  Center 
Space  Division,  Air  Force  Systems  Command 
Edwards  Air  Force  Base, 

California  93523-5000 


NOTICE  TO  ACCOMPANY  THE  DISSEMINATION  OF 
EXPORT-CONTROLLED  TECHNICAL  DATA 


1.  Export  of  information  contained  herein,  which  includes,  in  some  circumstances,  release  to 
foreign  nationals  within  the  United  States,  without  first  obtaining  approval  or  license  from  the 
Department  of  State  for  items  controlled  by  the  International  Traffic  in  Arms  Regulations  (ITAR), 
or  the  Department  of  Commerce  for  items  controlled  by  the  Export  Administration  Regulations  (EAR), 
may  constitute  a  violation  of  law. 

2.  Under  22  U.S.C.  2778  the  penalty  for  unlawful  export  of  items  or  information  controlled 
under  the  ITAR  is  up  to  2  years  imprisonment,  or  a  fine  of  $100,000,  or  both.  Under  SO  U.S.C., 
Appendix  2410,  the  penalty  for  unlawful  export  of  items  or  information  controlled  under  the  EAR 
is  a  fine  of  up  to  SI, 000,000,  or  five  times  the  value  of  the  exports,  whichever  is  greater;  or  for  an 
individual,  imprisonment  of  up  to  10  years,  or  a  fine  of  up  to  S2S0.000,  or  both. 

3.  In  accordance  with  your  certification  that  establishes  you  as  a  “qualified  U.S.  contractor,” 
unauthorized  dissemination  of  this  information  is  prohibited  and  may  result  in  disqualification  as 
a  qualified  U.S.  contractor,  and  may  be  considered  in  determining  your  eligibility  for  future  contracts 
with  the  Department  of  Defense. 

4.  The  U.S.  Government  assumes  no  liability  for  direct  patent  infringement,  or  contributory 
patent  infringement  or  misuse  of  technical  data. 

5.  The  U.  S.  Government  does  not  warrant  the  adequacy,  accuracy,  currency,  or  completeness 
of  the  technical  data. 

6.  The  U.S.  Government  assumes  no  liability  for  loss,  damage,  or  injury  resulting  from 
manufacture  or  use  for  any  purpose  of  any  product,  article,  system,  or  material  involving  reliance 
upon  any  or  all  technical  data  furnished  in  response  to  the  request  for  technical  data. 

7.  If  the  technical  data  furnished  by  the  Government  will  be  used  for  commercial  manufacturing 
or  other  profit  potential,  a  license  for  such  use  may  be  necessary.  Any  payments  made  in  support 
of  the  request  for  data  do  not  include  or  involve  any  license  rights. 

8.  A  copy  of  this  notice  shall  be  provided  with  any  partial  or  complete  reproduction  of  these 
data  that  are  provided  to  qualified  U.S.  contractors. 


FOREWORD 


This  finaj  report  on  Crack  Opening  and  Extension  in  Inert  Solid 
Propellant  was  submitted  by  Virginia  Polytechnic  Institute  &  State  University 
under  contract  number  F04611-85-C-0098  with  the  Air  Force  Astronautics 
Laboratory  (AFAL),  Edwards  Air  Force  Base,  CA  93523-5000.  W.  R.  Lloyd,  M. 
Rezvani,  T.  J.  Thelss,  Y.  Guo,  B.  Han,  and  J.  Jo  provided  valuable  assistance 
in  conducting  the  experiments  reported  in  this  study.  The  AFAL  project 
manager  was  Dr  C.  T.  Liu. 


This  technical  report  has  been  reviewed  and  Is  approved  for  publication 
and  distribution  In  accordance  with  the  distribution  statement  on  the  cover 
and  on  the  DD  Form  1473. 


Project  Manager  Chief,  Propellant  Development 

Branch 

FOR  THE  COMMANDER 


KENNETH  E.  MURPHY,  It  Col, 


Deputy  Director 
Solid  Rocket  Division 


j  Accession  For 

[  NTIS  GRA&I 
|  DTIC  TAB 
I  Unannounced 
!  Justlf ication_ 


□ 

□ 


% 


A 


By - - 

Distribution/ 
Avulia!  *  llty  r-.-i*. 
•Avail  nztC /  -  r 

'Diet  |  Special 


ICUniTV  CLASSIFICATION  OF  THIS  FAOC 


Block  11.  Experimental  Study 


Block  19. 


utilized  to  infer  the  distribution  of  the  dominant  eigenvalue  through  the  thickness  of  the 
model  in  a  cantinuun  sense.  Free  surface  values  compare  favorably  with  Ben  them's  analytical 
result.  The  near  tip  material  response  was  strongly  influenced  by  its  heterogeneous  nature. 
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;  INTRODUCTION 

X 

>When  cracks  occur  In  rocket  motor  grains,  they  may  significantly 
alter  the  subsequent  motor  performance.  Consequently  as  performance 
requirements  are  advanced.  It  becomes  critical  to  Incorporate  fracture 
control  concepts  Into  the  evaluation  of  the  structural  Integrity  of 
motor  grains.  While  most  such  concepts  are  typically  based  upon  two 
dimensional  analysis,  the  motor  geometry  and  conditions  of  use  suggest 
the  presence  of  significant  three  dimensional  effects.  Over  the  past 
decade,  attempts  have  been  made  to  develop  mathematical  models  for  pre¬ 
dicting  fracture  behavior  in  rocket  motors  but  model  validation  has  been 
Inhibited  by  lack  of  an  extensive  data  base  for  validating  three  dimen¬ 
sional  effects.  Moreover,  certain  such  effects  have  only  recently  been 
Identified,  and  have  not  been  evaluated  for  rocket  propellant. 

In  particular,  the  near  tip  behavior  of  cracks  In  solid  propellant 
material  during  opening  and  extension  of  the  crack  needs  to  be  quan¬ 
tified  experimentally.^ _ 

Program  Objectives 

The  objectives  of  the  present  program  were: 

I)  To  measure  the  components  of  the  displacement  field  on  the 
specimen  surfaces  around  the  crack  tip  at  time  Intervals  at  constant 
head  rate  extension  during  the  opening  and  extension  of  the  crack  In 
"biaxial ■  specimens  of  Inert  propellant  and  extract  surface  strain 
fields  from  the  data. 

II)  To  determine  the  dominant  singularity  order  for  the  crack  at 


the  specimen  surface  In  the  Inert  propellant  and  Infer  Its  distribution 
through  the  specimen  thickness  from  frozen  stress  measurements. 
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v)  Compare  the  the  thickness  averaged  results  with  numerical 
solutions  from  LEFM  and  surface  values  of  xu  with  analytical 
results. 

The  entire  test  procedure  described  above  was  first  applied  to  the 

large  specimens.  However,  since  the  primary  Interest  here  was  focused 

on  the  small  specimens,  we  shall  report  only  results  from  the  large 

specimens  and  then  present  details  In  conjunction  with  the  discussion  of 

•  • 
the  test  results  from  the  small  specimen  tests. 

From  the  through  thickness  photoelastic  tests,  stress  Intensity 
values  obtained  from  the  LEFM  method  of  analysis  described  In  Appendix  B 
agreed  closely  with  numerical  results  reported  In  Ref.  [51  for  large 
plates  with  through  cracks,  thus  confirming  that  xQ  »  |l-xu|  *  1/2  for 
the  through  thickness  data.  From  the  moire  analysis  at  global  strain 
levels  ranging  up  to  1%,  the  value  of  xu  (Eq.  B9  Appendix  B)  at  the  free 
surface  was  found  to  be  »  0.65  ±  .02,  agreeing  reasonably  well  with 
Benthem's  result  of  0.67.  Only  data  prior  to  crack  growth  were  used  In 
these  determinations.  Crack  growth  began  at  a  global  strain  level  of 
about  1.5%  after  which  xu  values  decreased. 

Typical  near  tip  through  thickness  photoelastic  and  free  surface 
moire  patterns  are  shown  in  Figs.  3a  and  3b  respectively  for  a  small 
specimen.  Analysis  of  digitized  photoelastic  data  from  several  global 
strain  levels  on  the  small  specimen  using  LEFM  from  which  the  thickness 
averaged  SIF  Is  obtained  using  the  data  zone  which  agrees  with  numerical 
LEFM  results  revealed  a  constricted  local  data  zone. 

Figure  4  Illustrates  how  xu  Is  determined  from  the  moire  data  using 
the  displacement  algorithm  described  In  Eq.  B9,  Appendix  B.  Tests  were 
run  over  a  global  strain  level  range  of  2.0  to  2.8%  without  crack  growth. 


The  moire  data  yielded  (xu)avg  *  0.65  which  again  agrees  well  with 
Bent hem's  result. 

It  Is  Important  to  note  that  the  data  zone  dominated  by  the  lowest 
eigenvalue  was  found  to  be  severely  constricted  In  the  small  specimen. 
For  example.  In  the  large  specimens,  the  moire  data  zone  extended  from 
0.60  to  3.5  mm  from  the  crack  tip  along  e  ■  */2.  However,  In  the  small 
specimen  with  constant  displacement  loading,  the  zone  was  constricted  to 
approximately  0.2  to  0.9  ran  from  the  crack  tip  (along  e  «  */2).  A 
similar,  though  less  severe,  contraction  occured  In  the  photoelastic 
data. 

Several  additional  tests  were  conducted  In  order  to  obtain  pre¬ 
liminary  findings  on  other  effects.  As  a  result  of  these  additional 
tests,  the  following  observations  are  made: 

1)  Increasing  model  thicknesses  from  6.4  mm  to  12.7  mm  did  not 
alter  values  of  xu  of  the  large  models. 

11)  Measurements  on  grown  cracks  suggested  the  presence  of  some 
crack  front  curvature  as  shown  In  Fig.  5. 

Ill)  Crack  tip  profiles  retained  a  nearly  elliptic  shape  as  shown 
In  Fig.  3b  and  revealed  no  macroscopic  evidence  of  crazing  In 
either  the  opened  crack  without  growth  or  grown  cracks. 

The  above  results  suggest  that: 

1)  The  zone  dominated  by  the  lowest  eigenvalue  Is  severely  con¬ 
stricted  In  the  small  biaxial  specimen. 


11)  The  Influence  of  the  free  surface  on  singularity  order  when 


using  thickness  averaged  data  Is  negligible. 

Ill)  In  nearly  Incompressible  materials  such  as  polyurethane,  the 
above  noted  effect  Is  quite  strong,  as  predicted  analytically 
by  Benthem  when  surface  measurements  are  made. 

Singularity  Order  Distribution 

In  addition  to  the  measurement  of  the  singularity  order  at  the  free 
surface.  It  was  desired  to  estimate  the  singularity  order  distribution 
through  the  specimen  thickness.  Having  previously  established  that  the 
body  geometry.  If  sufficiently  thick,  did  not  alter  the  value  of  the 
dominant  singularity  at  the  free  surface  [2],  one  might  expect  that  the 
distribution  of  the  dominant  eigenvalue  through  the  same  thickness  of 
body  might  also  be  unaffected  by  body  geometry  at  least  for  straight 
front  cracks.  Following  this  line  of  reasoning,  a  frozen  stress  photo¬ 
elastic  analysis  of  a  compact  bending  specimen  was  conducted  In  order  to 
estimate  the  distribution  of  the  dominant  stress  eigenvalue  \g  across 
the  specimen  thickness  (xq  ■  | xu-l | ) . 

A  compact  bending  model  consisting  of  the  geometry  pictured  In  Fig. 
6  was  cut  from  a  plate  of  stress  free  transparent  stress  freezing 
material.  A  straight  front  artificial  crack  was  machined  Into  the  model 
and  It  was  subjected  to  a  stress  freezing  cycle  under  four  point 
bending.  Subsequent  to  load  removal  at  room  temperature,  thin  slices 
mutually  orthogonal  to  the  crack  surface  and  the. crack  front  were 
removed  and  analyzed  utilizing  both  the  Tardy  and  fringe  multiplication 
analysis  techniques. 
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Section  A -A 


W  =  25.70  mm  *  =  50.8  mm 

d  =  6.35  mm  L  =  76.2  mm 

Q  s  12.82  mm  S  =  279.4  mm 

B  =  13.33  mm  P  =  2.32  N 

e  =  12.7  mm 


Raw  data  from  a  single  slice  are  shown  In  Fig.  7.  Algorithms  for 
converting  photoelastic  data  Into  fracture  parameters  are  discussed  In 
Appendices  B  and  C.  By  fitting  Eq.  B4  to  this  data,  the  slope  of  the 
linear  part  of  the  curve,  or  o°  can  be  determined.  This,  In  effect 
determines  (Eq.  C7)  for  the  case  where  xg  -  1/2  for  each  slice.  For 
slices  near  the  boundary,  computed  tq  values  begin  to  decrease. 

Beginning  at  this  location  from  the  free  surface,  we  arbitrarily  Impose 
a  linear  decrease  on  the  computed  value  of  down  to  zero  at  the  free 
surface.  By  Imposing  this  condition,  our  resulting  value  of  xq, 
measured  from  a  plot  of  Eq.  C9  agrees  with  both  Benthem's  results  and 
those  obtained  from  a  moire  Interferometric  analysis  at  the  free  surface 
[21.  A  log-log  plot  of  -r vs.  r  for  one  slice  near  the  center  of 
the  beam  Is  shown  In  Fig.  8. 

Figure  9  shows  the  values  of  across  the  thickness  of  a  compact 
bending  specimen  as  obtained  from  two  separate  tests  using  the  above 
procedure.  This  plot  reveals  the  thickness  of  the  boundary  layer 
through  which  xg  varies  for  an  Incompressible  homogeneous  material. 

Inert  Propellant  Tests:  Deformation  Measurement.  Methods  and  Results 

Specimen  and  Loading 

The  specimen  geometry  Is  Illustrated  In  Fig.  lb.  The  same  dimen¬ 
sions  were  used,  except  the  Inert  propellant  specimens  were  approxi¬ 
mately  0.2  In.  (5  mm)  and  0.6  In.  (15  mm)  thick.  The  artificial  cracks 
were  2a  ■  38  mm  long,  made  by  pressing  a  specially  prepared  razor  blade 
through  the  specimen  thickness;  the  blade  was  modified  by  grinding  Its 
short  edges  to  form  sharp  cutting  edges. 

The  specimen  material  Is  manufactured  by  Morton  Thlokol,  Inc., 
under  the  trade  name  H-24  Inert  Propellant.  The  principal  Ingredients 
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Fig.  9  Distribution  of  x  Across  the  Half  Thickness  of  a  Compact 
tending  Specimen.0 


are: 


polybutadiene  resin 

10%  by  weight 

potassium  chloride 

70% 

powdered  aluminum 

15% 

ammonium  perchlorate 

5% 

It  Is  a  dark  gray  rubbery  solid.  It  was  supplied  as  5  hid  and  IS  mm 
sheets,  presumably  cast  to  final  thickness  In  smooth  molds.  One  surface 
of  each  sheet  was  tacky  (or  sticky)  and  the  other  surface  was  not. 

The  material  is  a  heterogeneous  composite  of  very  stiff  Inclusions 
in  a  compliant  matrix  of  polybutadiene  rubber.  The  volume  fraction  of 
Inclusions  Is  approximately  75%.  The  principal  Inclusion,  potassium 
chloride.  Is  a  salt  with  essentially  cubical  shaped  grains  of  widths 
said  to  range  between  10  and  200  u*;  In  fact,  grains  as  large  as  500  ym 
(0.5  mm)  were  found,  and  a  great  abundance  of  400  ym  grains  were  found 
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In  the  material.  Many  crystals  of  potassium  chloride  were  broken  and 
they  exhibited  Irregular  shapes  Instead  of  cubes.  The  secondary 
Inclusions  were  small  compared  to  the  potassium  chloride. 

The  specimens  were  cemented  Into  relatively  rigid  aluminum  grips 
(Fig.  lb)  and  Installed  In  a  tensile  testing  machine.  The  Installation 
assured  parallel  translation  of  the  grips  (l.e.,  no  rotation  of  the 
grips).  Various  constant  cross-head  speeds,  or  grip  separation  rates, 
were  used  from  0.01  to  1.0  Inches  (0.25  to  25  mm)  per  minute.  A  typical 
graph  of  load  vs.  grip  extension  Is  shown  In  Fig.  10.  The  maximum  load 
was  reached  before  visible  crack  extension  occurred,  and  the  load 
continued  to  decrease  with  crack  growth.  Since  cross-head  speed  Is 
constant  for  each  test,  the  graph  also  represents  load  vs.  time. 

Displacement  Measurements,  Methods 

The  grid  method  was  used  for  In-plane  displacement  measurements  U 
and  V.  The  shadow  moire  method  was  used  simultaneously  on  the  back 
surface  of  the  specimen  to  measure  out-of -plane  displacements  W. 
Rationale  for  selecting  these  methods  In  preference  to  other 
experimental  techniques  Is  given  In  Appendix  D.  Details  of  the  grid 
method,  as  applied  here,  are  given  In  Appendix  E,  which  Includes  methods 
for  applying  the  grid  to  the  specimen.  Details  of  the  shadow  moire 
method  are  given  In  Appendix  F. 

A  separate  method  was  developed  and  demonstrated  for  acquiring  U,  V 
and  W  displacements  simultaneously  from  the  same  zone  of  observation, 
l.e.,  from  one  side  of  the  specimen.  This  Is  described  In  Appendix  G. 

It  utilizes  the  grid  method  for  U  and  V  and  two-beam  coherent  projection 
moire  for  W.  The  data  are  collected  by  a  single  photograph  from  which 


£0, Global  Strain 

2.5%  5.0%  7.5% 


the  grid  lines  can  be  digitized  and  the  W  displacement  field  extracted 
by  an  optical  filtering  process. 

Utility  of  W  Measurements 

The  algorithms  used  to  determine  fracture  parameters  In  this  work 
require  V,  alone,  and  not  W.  A  question  can  be  raised,  however, 
equivalent  to  that  for  assessing  the  magnitudes  of  finite  strains:  to 
what  extent  do  out-of-plane  rotations  Influence  the  change  of  length,  or 
relative  displacement  of  end  points,  of  a  line  segment?  If  W  Is 
Insignificant  In  the  calculation  of  In-plane  strains,  the  Influence  of  W 
on  the  fracture  algorithm  Is  negligible. 

That  Is  the  condition  realized  In  these  tests.  The  Influence  of  W 
on  ey,  for  example.  Is  given  by  the  Lagranglan  large  strain  relationship 

L  ,,  .  o  aV  .  r3Ui2  .  f 3V-j2  .  /-aW^  .  ,n 

'y  n  2  W  +  +  W  *  (iyJ  *  1  tl) 

where  L  signifies  the  Lagranglan  formulation.  Data  extracted  from  the 
W-dlsplacement  fields,  e.g.,  from  Fig.  11  or  12,  discussed  in  subsequent 
sections,  assures  that  (aW/ay)4,  Is  negligible  In  comparison  to  2aV/ay. 

This  Independence  of  W  was  observed  from  preliminary  tests. 
Accordingly,  the  utility  of  W  measurements  was  subordinated  to 
qualitative  observation  of  Its  distribution  during  the  crack  opening  and 
crack  growth  phases.  The  qualitative  character  could  be  assessed  from 
the  back  side  of  the  specimen  as  well  as  the  front.  This  removed  the 
requirement  for  observing  U,  V  and  W  at  the  same  location;  It  allowed  an 
easier  test  procedure  In  which  the  grid  method  and  shadow  moire  method 
could  be  set  up  Independently  of  each  other  on  front  and  back  surfaces, 
respectively. 


Fig.  11  Shadow  moire  pattern  of  W  field.  Test  8;  0.6  In _ 

thick  specimen;  1.0  In.  (25  wn)  per  win  loading  rate;  maximum  crac 
opening  prior  to  crack  growth;  sensitivity  O.OOZ  in.  (O.OS  mm)  per 


For  problems  where  the  Influence  of  W  Is  significant,  the  method 
described  In  Appendix  G  could  be  used. 

Test  Conditions 

Numerous  preliminary  tests  were  conducted  to  establish  the 
experimental  techniques.  Subsequently,  data  were  extracted  from  tests 
outlined  In  Table  2.  Test  No.  4  was  conducted  but  It  was  not 
successful. 

The  tests  were  conducted  at  room  temperature,  approximately  72°F. 
The  specimen  material  was  stored  In  Its  original  moisture-proof 
wrappings.  It  was  exposed  to  the  environment  while  the  specimen  was  cut 
to  shape,  while  the  grid  was  being  applied,  while  It  was  being  set  up 
for  testing  and  while  It  was  being  tested  --  In  all,  for  a  period  of  4 
to  8  hours.  Between  steps  In  this  process  It  was  protected  again  In 
moisture-proof  wrappings.  This  procedure  was  used  for  each  specimen 
except  No.  1,  which  was  exposed  to  the  environment  for  a  longer  time. 

Results:  Representative  Raw  Data 

A  sequence  of  photographs  of  the  crack  and  the  125  lines/inch  (5 
lines/mm)  grid  Is  given  In  Fig.  13  for  Test  No.  1.  Figure  13a 
Illustrates  the  crack  opening  phase,  In  which  the  width  of  the  crack 
Increases  as  the  grip  separation  proceeds.  Figure  13b  Illustrates  the 
crack  growth  phase.  For  each  photograph,  the  corresponding  load,  grip 
displacement  and  global  strain  Is  Indicated  In  Fig.  10  by  the  same  lower 
case  letter  as  used  In  Fig.  13.  Global  strain  eQ  Is  defined  as 
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where  6  Is  the  grip  displacement  and  hg  Is  the  Initial  height  of  the 
specimen. 

A  contour  map  of  out-of-plane  displacements  W  Is  Illustrated  by 
Fig.  11.  This  was  taken  In  Test  No.  8  with  a  0.6  Inch  (IS  mm)  thick 
specimen  and  1  Inch/mln  (25  mm/mln)  loading  rate.  It  shows  the 
condition  of  maximum  crack  opening  prior  to  crack  extension.  The 
contour  Interval,  or  sensitivity.  Is  0.002  Inches  (0.05  mm)  per  fringe 
order. 

Figure  12  Is  a  contour  map  of  W  obtained  by  the  simultaneous  U,V,W 
method,  using  a  0.6  Inch  (15  mm)  thick  specimen  In  the  crack  extension 
phase.  The  contour  Interval  Is  0.0012  Inches  (0.03  mm)  per  fringe 
order. 

Results:  Crack  Blunting 

Insight  Into  the  heterogeneous  nature  of  the  material  and 
probabilistic  character  of  Its  performance  can  be  found  In  Fig.  14  and 
Table  3.  Figure  14  shows  the  crack  and  the  surrounding  deformed  grid 
for  various  test  specimens  at  the  maximum  crack  opening  prior  to  crack 
growth.  The  degree  of  blunting  Is  rather  different  In  each  case, 
different  In  a  nonsystematlc  distribution.  The  cracks  terminate  with 
fairly  straight  blunted  tips,  essentially  perpendicular  to  the  crack 
length.  The  degree  of  blunting  will  be  quantified  here  by  the  length  of 
this  straight  portion,  called  blunt-tlp  height.  While  this  length  Is 
not  rigorously  defined  because  the  crack  tips  are  slightly  curved.  It 
can  be  measured  with  suitable  accuracy  as  the  distance  between  boundary 
points  where  the  radius  of  curvature  of  the  boundary  becomes  small. 

The  blunt-tlp  height  Is  given  In  Table  3.  The  three  replicate 
specimens  with  0.6  In.  (15  mm)  thickness  and  0.1  In. /min  (2.5  mm/mln) 
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Spec.  3. 


loading  rate  exhibit  significantly  different  blunting  behavior  both 
before  and  during  crack  growth.  It  suggests  that  the  local  arrangement 
of  Inclusions  at  the  crack  tip  plays  a  role  In  the  blunting  phenomenon. 

In  addition,  blunt-tlp  heights  measured  on  the  tacky  side  of  the 
specimen  regularly  exceeded  those  on  the  non-tacky  side.  The  average 
blunt-tlp  height  was  twice  as  great  on  the  tacky  side. 

Table  3  Indicates  a  typically  smaller  blunting  at  the  1.0  In. /min 
(25  mm/mln)  loading  rate.  This  Is  consistent  with  the  viscoelastic 
nature  of  the  rubber  matrix.  Since  these  results  were  all  from  the  non- 
tacky  side,  and  In  view  of  the  large  variation  found  with  replicate 
specimens,  the  results  cannot  be  viewed  as  a  definitive  trend. 

Crack-tip  blunting  that  occurs  during  crack  extension  Is 
Illustrated  in  Fig.  13b  and  also  by  several  examples  In  Figs.  15a  and 
b.  The  blunting  qualities  observed  during  the  crack  opening  phase  are 
exhibited  in  the  crack  growth  phase,  too.  The  magnitude  of  blunting 
varies  from  specimen  to  specimen.  It  Is  different  on  the  two  outside 
surfaces  of  the  specimen  at  any  Instant,  and  It  varies  Irregularly  with 
position  In  the  specimen  as  the  crack  front  advances.  Since  the  local 
arrangement  of  Inclusions  Is  a  random  function  of  position,  crack-tip 
blunting  might  be  expected  to  vary  with  position  of  the  advancing  crack 
tip.  With  regard  to  Its  magnitude,  the  results  Indicate  smaller  crack- 
tip  blunting  during  crack  growth,  compared  to  crack  opening,  when  the 
loading  rate  was  Increased  to  1.0  In. /min  (25  mm/mln). 

Results:  Fracture  Surfaces 

The  fracture  surfaces  were  very  rough.  Abrupt  changes  of  slope  of 
the  fracture  surface  were  evident,  both  along  the  crack  length  and  along 


irowth  phase.  Grid  500  lines  per 


the  crack  border,  l.e.,  across  the  specimen  thickness.  The  fracture 


appeared  to  be  confined  to  the  rubber  matrix.  Inclusions  were  not 


dislodged,  but  remained  Imbedded  In  the  matrix  on  both  sides  of  the 


fracture,  again  with  very  few  exceptions. 


Results:  Global  Deformations 


By  mutual  agreement,  displacement  and  strain  data  are  being 


transmitted  to  the  sponsor  on  computer  diskettes  for  the  tests  and  frame 


numbers  shown  In  Table  2.  The  results  of  data  reduction  are  given  for  a 


representative  case  In  the  sequence  of  graphs  of  Fig.  16a-h.  The  data 


are  taken  from  the  deformed  grid  of  Fig.  13a  [Image  (g)l,  which 


represents  the  maximum  crack  opening  prior  to  crack  growth. 


The  locations  of  grid  Intersection  points  were  digitized  and  their 


displacements  relative  to  the  undeformed  grid  [Fig.  13a,  Image  (a)l  were 


calculated.  Four-point  smoothing  was  used  (Appendix  E),  which  reduced 


the  experimental  error  of  reading  the  grid  and  It  attenuated  local 


anomalies  that  existed  In  the  real  deformation  of  the  heterogeneous 


specimen.  It  Is  a  smoothing  function  of  modest  strength,  which  does  not 


mask  the  nonuniform  deformations,  but  attenuates  them. 


Figure  16a  Indicates  the  displaced  positions  of  the  digitized  grid 


points.  Its  particular  purpose,  here.  Is  to  show  that  data  were  taken 


at  small  Intervals  near  the  crack  and  at  larger  Intervals  In  remote 


zones.  The  corresponding  locations  of  data  points  are  marked  In  the 


subsequent  graphs.  However,  they  are  shown  In  their  original  positions. 


not  their  displaced  positions.  The  x,y  positions  of  all  points  In  the 


subsequent  graphs,  and  the  displacements  and  strains  at  these  points. 


are  all  mapped  on  the  original  coordinates  of  the  specimen.  The  lips  of 


the  crack  are  represented  In  Fig.  16  by  a  single  line  of  zero  width. 


lacewent.  Test  1.  frame  12. — U  »  0.01_N 
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Test  1.  frame  12. 


Figure  16b  Is  the  V-dlsplacement  field,  obtained  by  connecting 
points  of  equal  V  displacement.  The  contour  Interval  Is  0.05  mm  or 
0.002  In.  and  the  displacement  at  any  point  Is  the  contour  interval 
multiplied  by  the  contour  number  N  at  the  point;  N  can  be  Interpolated 
between  curves.  Figure  16c  Is  a  magnified  view  of  the  V  field  near  the 
crack.  Some  lack  of  symmetry  and  lack  of  smoothness  Is  evident  In  both 
figures. 

What  Is  the  meaning  of  the  broad  black  region  near  the  crack?  It 
Is  not  the  shape  of  the  crack  opening.  Instead,  It  represents  the  zone 
around  the  crack  where  meaningful  data  are  not  available,  approximately 
one  grid  pitch  from  the  crack  lips. 

The  U  displacements  —  perpendicular  to  the  loading  direction  — 
are  plotted  in  Fig.  16d.  Note  that  the  contour  Interval  Is  very  small, 
1/5  that  of  the  V  contours.  Accordingly  the  gradients  of  U  are 
relatively  small  and  the  related  strain  contributions  are  small.  Note, 
too,  that  the  contours  are  free  of  severe  oscillations.  This  means  that 
the  effective  resolution  of  the  technique  Is  a  fraction  of  the  contour 
Interval,  or  about  2.5  ym  or  100  yin.  By  this  argument.  Irregularities 
In  a  contour  that  exceeds  about  1/3  the  distance  to  the  neighboring 
contour  represents  a  real  Irregularity  in  the  specimen  deformation;  in 
the  preceding  V  contours.  Irregularities  of  1/15  of  the  contour  spacing 
represent  real  behavior.  Curiously,  the  100  yin.  (2.5  ym)  resolution  Is 
equal  to  the  resolution  of  the  grid  digitizing  equipment  (Appendix  D), 
which  Implies  that  the  error  In  locating  grid  Intersection  points  did 
not  exceed  the  system  resolution. 

Figures  16e  and  f  are  contour  maps  of  the  dominant  strains,  cy. 
These  contours  are  much  less  smooth  than  the  displacement  contours. 


This  Is  natural,  because  strains  are  highly  sensitive  to  small  Irregu¬ 
larities  of  displacements.  It  Is  reasonable  to  expect  that  a  portion  of 
the  Irregularities  stems  from  experimental  error,  but  the  results  are 
reasonable  assessments  of  actual  specimen  behavior.  The  major  Irregu¬ 
larities  represent  local  strain  anomalies  in  this  heterogeneous 
material. 

These  strains  were  calculated  by  the  relationship  =  aV/ay,  which 

results  from  Eq.  1  when  the  derivatives  of  U  and  W  are  negligible. 

« 

The  calculation  of  strains  did  not  use  the  smoothed  displacement 

data.  Instead,  for  each  array  of  4  measurement  points  forming  a 

rectangle,  the  change  of  length  of  the  vertical  sides  of  the  rectangle 

were  determined  and  averaged;  similarly  for  the  horizontal  sides;  the 

three  strains,  e  ,  e  and  y  t  were  determined  from  these  averaged 
y  x  xy 

changes  of  length  and  assigned  to  the  center-point  of  the  rectangle. 

This  Interpretation  of  the  raw  data  afforded  only  a  small  degree  of 
smoothing,  providing  some  attenuation  of  severe  irregularities.  Another 
factor  that  provided  attenuation  of  irregularities  was  Increase  of  gage 
length.  The  array  of  4  data  points  forming  the  rectangle  was  not 
necessarily  the  4  nearest  points;  the  rectangles  used  for  strain 
calculations  were  larger,  so  the  effective  gage  lengths  —  the  sides  of 
the  rectangle  —  were  larger  and  the  local  Irregularities  were 
smoothed.  The  locations  of  the  data  points  used  for  the  calculations 
are  printed  on  each  figure.  In  these  experiments,  the  pitch  of- the  grid 
printed  on  the  specimen  was  less  than  half  the  width  of  the  largest 
(rigid)  Inclusions,  and  extensive  small-scale  Irregularities  of 
deformation  are  Inevitable.  The  data  reduction  techniques  were  designed 
to  mask  or  attenuate  some  of  these,  for  otherwise  the  contour  maps  of 


strain  would  show  great  oscillations  Instead  of  the  general  strain 
trends. 


Figure  16g  shows  the  contours  of  ex,  but  with  a  smaller  contour 
Interval.  Clearly  these  strains  are  very  small  and  Insignificant 
compared  to  ey. 

Figure  16h  shows  y  ,  the  In-plane  shear  strains  near  the  crack, 
xy 

The  shear  strains  are  significant  near  the  tip  of  the  crack,  and 
Insignificant  elsewhere.  This  is  also  evident  In  Fig.  13a  (Image  g) 
which  shows  the  change  of  grid  elements  from  rectangular  to  rhomboldal 
shape  near  the  crack  tip. 

The  shear  strain  was  calculated  by  the  strain  relationship 

Yxy  *  aU/ay  +  3V/3X  (3) 


A  more  rigorous  calculation  would  use  the  Lagranglan  large  strain 
relationship 


aU  +  aV 

arcs  In 


(1  +  ex)(l  +  cy) 


(4) 


For  shear  strains  near  0.3,  the  arcsin  and  the  angle  in  radians  are 

nearly  equal.  As  already  Indicated,  ex,  or  aU/ax,  Is  negligible  and  the 

product  of  derivatives  of  W  are  negligible.  In  the  region  where  y  Is 

xy 

large,  aU/ay  Is  much  smaller  than  aV/ax,  as  seen  In  the  grid  of  Fig.  13a 
(Image  g).  Under  these  conditions,  the  right  side  of  Eq.  4  reduces  to 
aV/ax,  so  use  of  Eq.  3  Is  justified  as  a  good  approximation. 

Figure  17a  shows  the  V-dlsplacement  fields  around  cracks  that  have 
reached  their  maximum  opening  prior  to  crack  growth.  The  vertical  line 
Indicates  the  location  of  the  crack  tip  In  each  case.  The  bar  In  each 


lower-left  corner  Is  a  scale  representing  i' mm  length  In  the  specimen. 
The  contour  Interval  is  K,  given  in  mm  of  displacement  per  contour 
level.  Again,  data  points  are  plotted  In  their  Initial  positions.  The 
contour  maps  show  essentially  the  same  global  features  for  each  test. 
Local  Irregularities  appear  in  the  contours  and  they  appear  to  be 
randomly  distributed. 

The  corresponding  strain  fields,  ey,  are  shown  in  Fig.  17b.  The 

crack  positions  and  scales  are  defined  as  before.  The  contour  Interval 

♦ 

K  Is  the  strain  (m/m)  per  contour  level.  Again  the  general  features  are 
the  same  for  each  test.  Local  random  Irregularities  are  extensive,  as 
expected  for  the  heterogeneous  specimen  material.  The  patterns  show 
(Irregular)  cusp-shaped  contours  symmetrical  with  the  crack  plane;  on 
lines  x  =  constant,  ahead  of  the  crack,  ey  reaches  Its  maximum  value  at 
some  distance  from  the  crack  plane. 

Strain  rate  effects  are  shown  In  Fig.  18.  The  rate  of  deformation 
Is  not  constant.  Instead,  as  the  specimen  Is  stretched  with  a  uniform 
rate  of  grip  displacement,  the  deformation  Increases  slowly  at  first  and 
then  more  quickly. 

In  the  three  time  Intervals  of  one  second  each,  the  crack  opening 
Increased  progressively  by  0.08,  0.16  and  0.22  mm.  This  strong 
nonlinearity  Indicates  a  redistribution  of  material  compliance  during 
the  test.  It  changes  from  a  uniform  compliance  throughout  the  specimen 
at  the  beginning  of  the  test  to  a  subsequent  condition  of  Increased 
compliance  locally  where  strains  become  high.  A  physical  model  of  such 
material  behavior  Is  proposed  In  Appendix  H. 


Interpretation:  Process  Zone 

The  model  proposed  In  Appendix  H  Is  developed  to  explain  the 

physical  behavior  In  the  process  zone,  or  pre-crack  zone.  The  model  Is 

explained  In  terms  of  Idealized  arrangements  of  rigid  Inclusions  In  the 

rubber  matrix,  but  the  general  features  should  apply  to  the  random 

arrangement  present  In  the  real  specimen  material.  Since  the  specimen 

has  random  qualities,  however,  the  model  represents  generalizations  or 

first  approximations  of  its  performance.  The  model  prescribes  the 

« 

following  performance. 

1.  The  process  or  pre-crack  zone  extends  as  a  triangular  region  ahead 
of  the  blunt  crack  tip.  Virtually  no  Inclusions  lie  In  this  zone. 
Instead,  the  zone  separates  Inclusions  that  were  previously  located  on 
each  side  of  a  very  narrow  path  of  matrix  rubber. 

2.  The  material  In  the  pre-crack  zone  Is  matrix  rubber  that  spans  the 
zone  as  narrow  ligaments.  The  magnitude  of  the  normal  strain  in  the 
ligaments  Is  of  the  order  of  *  5. 

3.  Let  the  rubber  surrounding  the  Inclusions  In  the  undeformed  state  be 
categorized  as  having  basically  vertical  elements  (with  their  long 
dimension  parallel  to  the  load  direction)  and  horizontal  elements.  The 
ligaments  are  extended  bodies  of  matrix  rubber  that  were  previously 
vertical  elements  between  inclusions. 

4.  As  the  deformation  develops  from  zero,  horizontal  elements  of  matrix 
rubber  sustain  the  stress.  The  corresponding  material  stiffness  Is 
high.  As  deformation  Increases,  voids  are  developed  between  horizontal 
surfaces  of  Inclusions  and  rubber.  Then  debonding  and  relative  sliding 
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between  rubber  and  inclusions  begins  at  the  vertical  Interfaces. 
Debonding  and  sliding  ceases  where  attached  horizontal  elements  of 
matrix  key  the  rubber  to  the  Inclusions.  Substantially  Increased 
compliance  accompanies  this  voiding  and  stretching  of  vertical  elements 
of  rubber. 

5.  The  stretched  vertical  elements  of  rubber  carry  the  load  In  both  the 
pre-crack  zone  and  the  surrounding  finite  strain  zone.  This  material 
stretches  to  form  ligaments  In  the  pre-crack  zone,  which  subsequently 
rupture.  Crack  growth  takes  place  as  the  ligaments  In  the  pre-crack 
zone  rupture. 

6.  Rubber-rich  surface  layers  bound  the  outside  surfaces  of  the 
triangular  pre-crack  zone.  They  stretch  to  span  the  pre-crack  zone. 

The  vertical  leg  of  the  triangular  surface  layer  Is  the  apparent  blunt 
end  of  the  crack. 

7.  Large  tensile  strains  and  shear  strains  are  developed  In  the  surface 
layer.  The  principal  tensile  strain  remains  nearly  vertical.  Small 
cracks  are  formed  ahead  of  the  main  crack,  presumably  wherever  the 
principal  strain  reaches  a  critical  value.  The  crack  planes  would  be 
essentially  horizontal  like  the  main  crack,  but  the  crack  lips  deform 
Into  tall  elliptical  shapes  In  the  severe  strain  field. 

8.  As  the  specimen  Is  stretched,  elliptical  voids  coalesce  and  also  the 
blunt  end  of  the  surface  layer  progressively  ruptures.  Thus,  the  crack 
tip  advances  as  new  pre-crack  material  Is  developed  ahead  of  It. 

9.  The  end  of  the  pre-crack  zone  In  the  surface  layer  probably  lags 
that  In  the  Interior. 


10.  At- any  junction,  the  preferred  path  of  the  crack  passes  through  the 
thinner  layer  of  matrix  rubber.  Inclusion  arrangements  that  Inhibit 
high  strains  can  act  as  crack  arresters.  Then  a  new  pre-crack  zone 
develops  at  a  neighboring  strain  concentration. 

Determination  of  xu  from  Cracked  Inert  Propellant 

Since  tests  on  large  polyurethane  plates  revealed  no  thickness 
effect  upon  the  determination  of  xu  for  6.4  and  12.7  mm  thick  specimens. 
It  was  decided  to  utilize  the  thin  (5  mm  thick)  plates  of  Inert 
propellant  mainly  to  make  specimens  for  pilot  tests  to  determine 
measurement  techniques  and  then  utilize  the  thicker  (15.25  mm)  plates 
mainly  for  determination  of  xu  values  at  various  displacement  head 
rates.  Using  this  approach,  values  of  xu  were  determined  at  several 
different  global  strain  levels  up  to  Initiation  of  growth  for  three 
different  head  displacement  rates. 

A  brief  description  of  the  data  collection  and  reduction  techniques 
used  will  now  be  presented.  Figure  15b  shows  the  surface  profile  of  an 
opened  crack  with  a  surrounding  grid  of  *  20  i/mm.  Data  are  taken  only 
along  the  vertical  lines  passing  through  the  crack  tip.  Data  points 
(l.e.  locus  of  the  Intersection  points  of  the  grating  lines)  are  col¬ 
lected  and  stored  In  a  microcomputer  using  a  digitizing  table  with  reso¬ 
lution  of  0.025  mm.  Appropriate  software  were  written  for  use  In  data 
collection  and  reduction.  The  error  Introduced  Into  the  data  due  to 
resolution  error  was  found  to  be  negligible  so  long  as  the  photographic 
prints  were  magnified  at  least  10X.  Data  were  collected  both  above  and 
below  the  crack  tip  and  were  then  averaged. 

The  coordinates  of  the  grid  points  In  the  no  load  condition  are 
necessary  for  data  reduction.  The  no  load  coordinates  can  be  found  in 


one  of  two  ways.  First,  the  enlarged  photograph  of  the  no-load  grating 
can  be  read  as  previously  described.  Alternatively,  the  no  load  data 
point  locations  can  be  computed  assuming  perfect  grid  application.  Both 
methods  were  used  and  gave  essentially  the  same  results.  The  present 
analysis  was  based  upon  the  perfect  grid  assumption. 

As  noted  In  the  previous  section,  the  cracks  In  the  Inert  propel¬ 
lant  biaxial  specimen  tests  exhibited  severe  blunting  In  some  cases  but 
not  In  others.  In  order  to  account  for  the  blunting  effect.  It  was 
necessary  to  further  generalize  the  algorithm  given  by  Eq.  B9  for  con¬ 
verting  displacement  data  Into  the  appropriate  fracture  parameter  (in 
this  case  xu).  Basically,  this  Involved  the  Inclusion  of  a  term  (V)0  In 
the  equation  for  the  displacement  V  along  e  =  ir/2.  Thus  Eq.  B8,  when  so 
modified  becomes: 

v  -  <V)o  +  DyrX“  (5) 

whence 

In  [V-(V)Q]  =  In  Dy  +  xyln  r  (6) 

The  value  of  V  Is  found  by  subtracting  the  y  coordinate  under  load 
from  Its  no-load  value.  Figure  19  Indicates  schematically  the  role  of 
(V)q  and  Figure  20  shows  test  data  from  which  xu  Is  determined.  It  Is 
Important  to  collect  data  over  a  much  larger  range  than  the  zone 
dominated  by  xy  so  as  to  clearly  delineate  this  data  zone. 

Oata  from  four  different  tests  were  used  to  determine  xu  values. 
Information  pertinent  to  these  tests  Is  presented  In  Table  4.  The  data 
from  the  specimens  utilizing  the  5  lines/mm  grating  were  obtained  by 
reading  both  the  top  and  bottom  of  the  horizontal  grid  lines  Instead  of 
the  central  values.  The  location  of  the  data  zone  dominated  by  the 
lowest  eigenvalue  xy  Is  seen  to  be  further  from  the  crack  tip  due  to 
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Table  4:  Test  Data  for  x  Determination 


Head  Rate  (iran/min)  Grating  Frequency  a/mm)  Data  Zone  (mm) 

.25  20  0.5-1. 5 
2.5  5  0.6-1. 6 
2.5  20  0. 6-1.0 
25  5  1. 0-2.0 


Table  5:  Test  Data  for  x  Measurements 


-  Range  (0.5-1.5]mm  Dlspl.  Rate  *  0.01  In/mln  or  0.25  In/min 


«(*) 

1.5 

0.625 

3.0 

0.687 

5.0 

0.778 

5.5 

0.748 

XUAVG  =  0.71 


-  Range  [0.6-1.6]mm  '  Dlspl.  Rate  =  0.1  In/min  or  2.5  mm/min 


1.25  0.689 
2.10  0.698 
2.92  0.688 


xyAVG  =0.69 


-  Range  [0.6-1.0]mm 


Dlspl.  Rate  =  0.1  In/min  or  2.5  mm/min 


z(%) 

\l 

1.25 

0.707 

1.67 

0.758 

2.10 

0,630 

2.92 

0.652 

XyAVG  =0.69 


-  Range  [1.0-2.0jmm 


Dlspl.  Rate  *  1.0  In/min  or  25  mm/min 


blunting  than  In  the  polyurethane  models  and  Its  location  varies 
somewhat  In  the  different  tests.  The  latter  is  attributed  mainly  to  the 
heterogeneity  of  the  Inert  propellant  although  the  displacement  rate  may 
have  had  some  Influence. 

The  resulting  values  of  x are  presented  in  Table  5.  The  average 
value  of  xu  appears  to  decrease  slightly  with  Increasing  global  dis¬ 
placement  rate.  However,  It  should  also  be  noted  that  some  crack  growth 
appeared  to  occur  for  the  extreme  load  rates  (i.e.  Tests  8  and  9). 
Consequently,  at  this  juncture,  it  would  appear  that  xu  Is  relatively 
Independent  of  the  head  displacement  rate.  Taken  collectively,  for  the 
15  measurements,  the  results  yield  an  average  value  of  xu  =  0.69  with  a 
standard  deviation  of  ±  0.06.  These  results  compare  favorably  with 
Benthem's  result  of  0.67,  suggesting  that  the  continuum  description  of 
this  phenomenon  Is  still  valid  during  the  crack  opening  process  in  the 
Inert  propellant. 

Returning  now  to  Figure  9,  and  plotting  x^  =  (1-xJ  =  0.31  avg,  one 
would  have  a  data  point  which  would  lie  approximately  on  the  lower  bound 
of  the  frozen  stress  test  data.  Since  the  thickness  of  the  compact 
bending  specimens  were  approximately  the  same  as  the  thicker  biaxial 
inert  propellant  specimens,  one  may  infer  the  transition  zone  for 
x  (or  x,,)  to  be  the  same  thickness  also.  Finally,  we  arrive  at  the  x 
distribution  for  the  Inert  propellant  by  adjusting  the  transition  zone 
to  agree  with  the  surface  value  of  0.31  and  this  is  shown  In  Fig.  21. 

In  view  of  the  heterogeneous  nature  of  the  inert  propellant.  It  appears 
Inappropriate  to  attempt  further  refinement  of  the  distribution. 

However,  It  should  be  noted  that  some  crack  growth  did  occur  in  the 
extremely  slow  and  fast  displacement  rate  cases,  producing  more  scatter 
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Fig.  21  \  Distribution  for  Inert  Propellant 


In  xu  for  those  tests  (8  and  9)  than  might  otherwise  have  been 
experienced.  In  any  event.  It  appears  that  the  values  of  \g  through  the 
transition  zone  In  the  Inert  propellant  biaxial  specimens  may  be 
slightly  lower  than  In  the  compact  bending  polyurethane  specimens. 

SUMMARY  AND  RECOMMENDATIONS 
Summary 

The  foregoing  study  consisted  of  four  main  parts.  These  parts  and 
their  findings  are  as  follows: 

I  -  Preliminary  tests  were  conducted  on  transparent  polyurethane 
specimens  to  evaluate  the  Influence  of  the  free  surface  effect  on  the 
order  of  the  singularity  of  the  biaxial  test  configuration  and  the  zone 
dominated  by  this  effect  together  with  the  Influence  of  specimen 
thickness  on  this  effect.  Results  showed  that 

I)  The  zone  dominated  by  the  free  surface  singularity  effect 
exists  but  Is  severely  reduced  by  the  biaxial  specimen  geometry  and 
extends  from  approximately  0.2  to  0.9  mm.  In  Isotropic,  homogeneous, 
nearly  incompressible  polyurethane. 

II)  The  Influence  of  the  free  surface  on  the  stress  singularity 
order  at  the  crack  tip  Is  negligible  for  through  thickness  data. 

III)  Additional  observations  revealed  that: 

a)  Increasing  specimen  thickness  from  6.4  to  12.7  mm  did 
not  alter  the  values  of  xu  In  large  polyurethane  models. 

b)  Some  crack  front  curvature  was  observed  when  crack 
growth  occurred. 

c)  Crack  tip  profiles  retained  a  nearly  semi-el llptlc  shape 


when  loaded  and  revealed  no  macroscopic  evidence  of  crazing  during 
either  opening  or  growth  stages. 


II  -  Preliminary  tests  were  conducted  on  cracked  biaxial  Inert 
propellant  specimens  In  order  to  determine  appropriate  methods  for 
measuring  near  tip  displacement  and  strain  fields  during  various  con¬ 
stant  cross-head  displacement  rates.  The  grid  method  was  selected  for 
in-plane  displacement  measurements  and  the  shadow  moire  method  for  out- 
of-plane  measurements.  Methods  for  printing  gratings  on  the  Inert 
propellent  surface,  and  Instrumentation  to  record  sequential  Images  of 
the  deformed  grid  and  shadow  moire  fringes  were  developed.  Techniques 
for  data  extraction  and  data  reduction  were  Implemented.  In  addition,  a 
test  method  for  simultaneous  determination  of  U,V  and  W  displacements 
from  a  single  photograph  was  demonstrated. 

III  -  Local  displacement  and  strain  fields  In  cracked  biaxial  inert 
propellant  specimens  were  determined  for  0.2  and  0.6  In.  (5  and  15  mm) 
thick  specimens  loaded  at  various  cross-head  rates  from  0.01  to  1.0  In. 
(0.25  to  25  mm)  per  minute.  Ho  systematic  dependence  on  thickness  or 
loading  rate  was  found.  The  character  of  the  global  fields  was 
established.  Local  Irregularities  associated  with  the  heterogeneous 
nature  of  the  specimen  material  were  abundant. 

Strain-rate  effects  were  Investigated  for  the  crack  opening  phase. 
The  deformation  accumulates  slowly  at  first  and  then  progressively  more 
rapidly  as  the  specimen  Is  stretched  at  a  uniform  rate.  This  nonlinear 
response  was  attributed  to  a  redistribution  of  local  material  compliance 
around  the  crack,  consistent  with  the  model  of  Appendix  H. 

IV  -  Results  of  extraction  of  free  surface  singularity  order  from 
measurements  on  cracked  biaxial  Inert  propellant  specimen  surfaces  and 
inference  of  Its  distribution  through  the  thickness  from  frozen  stress 
results  yielded: 


I)  An  average  value  of  xu  =  0.69  over  15  tests  with  a  standard 
deviation  of  0.06  which  Included  global  strain  levels  of  1.5  to  10%  and 
head  rates  of  0.25,  2.5  and  25  mm/min. 

II)  Using  frozen  stress  data  from  a  compact  bending  specimen  of 

approximately  the  same  thickness  as  the  biaxial  specimens,  a  transition 

zone  for  xg  was  found  to  be  about  20%  of  the  specimen  thickness  over 

which  x  varied  from  0.31  to  0.50. 
o 

These  results  reveal  that,  using  a  continuum  approach,  the  LEFM 
stress  singularity  is  lost  near  the  free  surface  of  the  crack  front/free 
surface  Intersection  point  but  Increases  through  a  transition  zone  to 
Its  LEFM  value. 

Since  the  value  of  the  lowest  dominant  eigenvalue  differs  from  that 
associated  with  classical  fracture  mechanics,  one  Is  led  to  address  the 
issue  of  how  the  quantified  knowledge  of  this  phenomena  can  be  Intro¬ 
duced  Into  the  framework  of  classical  fracture  mechanics  so  It  can  be 
utilized  In  engineering  analysis.  A  simple  way  In  which  to  achieve  this 
Is  to  introduce  a  stress  Intensity  factor  which  would  be  required  to 

produce  the  same  value  of  as  one  measures  In  the  actual  tests  where 

max 

\g  Is  not  equal  to  the  classical  value.  Calling  this  the  corresponding 
SIF,  one  may  write 


“oAd  (\> 


(7) 


Now,  substituting  relation  (7)  into  Eq.  (C6),  one  has: 


Tmax 


*g^cor^Ap 
/2tT  r^2 


(8) 


Now,  combining  Eqs.  (B3)  and  (8),  one  has 


whence 


Xo^Kcor^Ap  _ _ S<Ap 

✓27  r1^2  /87  r1^2 


^KcorUp 


(9) 


Since  values  of  K  and  Kcor  are  obtained  by  extrapolating  plots  of  KAp 
and  (Kcor) Ap  vs  across  a  nonlinear  zone  to  the  origin.  It  follows 
that. 


(10) 


It  should  be  emphasized  that  the  Introduction  of  relation  (7)  Is  quite 
arbitrary  and  Is  based  upon  purely  dimensional  grounds.  Moreover,  since 
the  analysis  here  Is  quasi-two  dimensional,  one  depends  on  moire  results 
using  Benthem's  analysis  to  validate  the  stress  algorithm  at  the  free 
surface. 

Equation  (9)  then,  may  be  applied  In  the  transition  zone  to  convert 
the  SIF  obtained  from  LEFM  Into  the  value  which  would  exist  there  had 
the  singularity  order  been  0.5.  In  the  present  case: 


K=or  -  JOT  '  !•«  (ID 

and  so  Kcor  may  be  viewed  as  a  very  conservative  way  of  treating  free 
surface  effects.  On  the  other  hand,  If  one  Is  dealing  with  thickness 
averaged  values  Instead  of  surface  values,  no  correction  for  the  effect 
would  appear  to  be  Indicated.  It  should  be  emphasized  that  this  discus¬ 
sion  Is  based  upon  a  continuum  approach. 


Recommendations 

The  foregoing  studies  revealed  that  the  location  of  the  zone  domin¬ 
ated  by  the  lowest  eigenvalue  was  further  from  the  crack  tip  in  Inert 
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propellant  than  In  polyurethane  and  that  significant  changes  of  material 
properties  occurred  in  the  process  zone  for  propellant  but  not  for  the 
polyurethane.  Moreover,  severe  crack  tip  blunting  occurred  In  the 
propellant.  Special  blunting  characteristics  were  associated  with  a 
rubber-rich  surface  layer  on  the  propellant  material.  Finally,  tests  on 
the  polyurethane  showed  some  thumbnalllng  when  the  cracks  grew.  It  is 
recommended  that  an  exploratory  study  be  undertaken  using  Inert  propel¬ 
lant  material,  binder  material  only  and  binder  material  with  additive  to 
determine  the  differences  between  crack  opening  and  growth  mechanisms  in 
these  materials  and  to  attempt  to  extract  values  from  such  tests 
employing  biaxial  specimens.  The  inert  propellant  material  should  be 
machined  from  bulk  stock,  rather  than  cast  In  plate  form,  to  eliminate 
the  effects  of  rubber-rich  surface  layers.  It  is  also  recommended  that 
a  frozen  stress  analysis  of  a  thumbnalled  crack  front  simulating  that 
exhibited  by  a  growing  crack  in  Inert  propellant  be  undertaken  and 
results  be  compared  with  straight  front  crack  results. 
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APPENDIX  A  INSERTION  OF  CRACKS  IN  POLYURETHANE  MODELS 


A  cutter  (Fig.  A-l)  made  by  shaping  and  sharpening  the  edges  of  a 
putty  knife  was  forced  through  the  thickness  of  the  model  material  to 
make  a  straight  front  crack.  Slotted  aluminum  blocks  were  used  above 
and  below  the  polyurethane  to  Insure  that  the  crack  fronts  were  normal 
to  the  specimen  surfaces. 


APPENDIX  B  ALGORITHMS  FOR  CONVERTNG  OPTICAL  DATA  INTO  FRACTURE 


PARAMETERS  (LEFM  AND  X..  UNDEFINED 


In  developing  algorithms  for  converting  free  surface  data  into 
fracture  parameters,  the  authors  used  both  linear  elastic  fracture  mech¬ 
anics  (LEFM)  concepts  and  the  results  of  Benthem's  analysis  (loc,  clt.). 
For  benefit  of  the  reader  and  preservation  of  continuity  both  will  be 
briefly  reviewed  here. 


Mode  I  LEFM  Algorithm  for  Photoelastic  Data 


We  may  describe  the  near  tip  stresses  In  LEFM  for  two  dimensional 
problems  as  (Ref.  [ B1 ] ) 


(2»r) 


C/2  ^Ij^^  ~ 


where  o^j  are  the  near  tip  In-plane  stress  components,  Kj  Is  the  Mode  I 
stress  Intensity  factor  (SIF),  o^j  may  be  regarded  as  the  near  tip  con¬ 
tribution  of  the  non-singular  stresses  which,  experience  shows,  can  be 
represented  as  a  set  of  constant  values  very  near  the  crack  tip  and 
(r,e)  are  polar  coordinates  as  shown  In  Fig.  B-l.  Along  e  =  w/2,  we  may 
compute  from  Eqs.  (Bl)  [when  (o°)2  Is  small  relative  to  8t£  [Bl]] 
where  o°//5  represents  the  contribution  of  the  non-singular  stress  state 


1 


to  locally  as: 
max 


max 


(8*r) 
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(82) 


Now  defining  an  "apparent"  SIF, 

<ki>ap  ■  W8*r>1/2  <B3> 

—  1/2  — 

and  normalizing  with  respect  to  o(*a)  '  where  a  represents  the  remote 
stress  and  a  the  half  crack  length,  we  have 


(K1^AP  s  K1  c!  /rxl/2 

o(*a)1/2  7(na)1/2  o  a 


(B4) 


which  suggests  a  linear  zone  In  a  plot  of 


^ki'ap 


vs  <*) 


1/2 


0 

with  a  slope  of 

o 


Experience  shows  this  zone  to  lie  between  (r/a)^2  values  of  approxi¬ 
mately  0.2  to  0.4  (or  above)  In  most  two  dimensional  problems.  By 
extracting  optical  data  from  this  zone  and  extrapolating  across  a  near 
tip  non-linear  zone,  an  accurate  estimate  of  K^/o(na)^2  can  be 
obtained.  This  Is  Illustrated  for  a  large  specimen  In  Fig.  B-2. 
Details  of  this  approach  are  found  In  Ref.  [B1J. 


Mode  I  LEFM  Moire  Algorithm 

Corresponding  to  Eqs.  (Bl)  we  have,  for  displacements  for  the  two 


dimensional  case  (Ref.  [B2]) 


o 


o 


o 


F19*  B2  Determination  of  Normalized  from  Photoelastic  Data 
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a 


U  *  CK1Gx(e)r1/2  -  f-  rcosa 
V  =  CKXG  (e)r1/2  -  rslne 


where  C  Is  a  constant  containing  the  elastic  constants,  o°  is  a  non¬ 
singular  stress  term  and  Es  is  the  shear  modulus. 

In  the  photoelastic  algorithm,  we  included  the  effect  of  the  non¬ 
singular  stress  since  It  was  Independent  of  r.  However,  In  Eqs.  (B5), 
we  may,  for  sufficiently  small  values  of  r,  neglect  the  effect  of  the 
a0  terms.  Thus,  along  e  =  w/2,  we  may  write. 


V  =  C'KjT  '  (B6) 

where  a  plot  of  V/r1/2  vs  r1/2  should  yield  a  horizontal  line  with 
ordinate  proportional  to  Kj.  Sometimes  however,  the  o°  term  does  make  a 
contribution  in  the  data  zone  so  we  must  use  the  form 

V  .  C"(Kj)Ap  r1/2  (B7) 

and  extrapolate  to  the  origin  as  before  a  plot  of  V/r^2  vs  r^2. 


Algorithms  for  Accounting  for  Free  Boundary  Effect  on  the  Crack  Tip 
Singularity  at  the  Free  Surface 

Following  the  reasoning  In  LEFM  and  focusing  only  on  the  leading 
term  In  Benthem's  series  expression  for  V,  we  may  write  directly, 
along  0  =  */2 


V  *  Dyr 


whence 


log  V  =  log  Dy  +  xjogr 


so  that  the  slope  in  the  linear  zone  of  a  plot  of  log  V  vs.  log  r  will 


yield  xu. 
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APPENDIX  C  ALGORITHMS  FOR  CONVERTING  OPTICAL  DATA  INTO  FRACTURE 
PARAMETERS  (x  UNDEFINED^ 

In  developing  algorithms  for  converting  frozen  stress  data  into 
values  of  x  ,  the  authors  used  both  the  linear  elastic  fracture  mech- 

o 

anics  (LEFM)  concepts  described  in  Appendix  B  and  the  results  of 
Benthem's  analysis  (loc,  ci t . ) .  The  developments  from  the  latter 
results  are  described  below. 


Photoelastic  Algorithms  for  Accounting  for  Free  Boundary  Effect  on 


the  Crack  Tip  Singularity  at  the  Free  Surface 


Although  the  above  noted  effect  involves  a  three  dimensional  state 
of  stress,  one  may  consider  measurements  on  the  free  surface  itself  and 
within  thin  layers  parallel  to  the  free  surface  as  being  associated  with 
a  quasi-two  dimensional  stress  state.  Thus,  one  may  begin,  as  In  LEFM, 
with  a  two  dimensional  stress  function  f (c)  the  form: 

Z  *  ^  (Cl) 


where  t  =  re  w  and  is  measured  from  an  origin  at  the  crack  tip  (Fig. 


Bl).  Expanding  f(c)  in  a  Maclaurin  series  about  the  crack  tip  and 


letting  ;  ♦  0  we  recover 


/2u  r  ° 

where  Kx  may  be  designated  as  a  stress  eigenfactor.  Now  substituting 
(C2)  into 


oxx  =  ReZ  -  ylmZ 


®yy  =  ReZ  +  ylmZ 


oxy--yReZ'. 


we  obtain,  using  o°  to  account  for  non-biaxiality  of  remote  loads. 


/2n  r 


—  {cosx  e  -  x  sine  sin(x  +  l)e}  -  ol 

X  a  a  a 

a 


^1 


Oyy  * - -  {cosxoe  +  X0s1ne  sin(xo  +  l)e> 

n*  r  ° 


o  * - -  {x  sine  cos(x  +  l)e> 

G  G 

/2n  r  ° 


where  o°  plays  the  same  role  as  In  LEFM.  From  Eqs.  (C3)  one  may 
compute 


3 


,,°yy  "  °xxx2  .  2  ,1/2 

W  =  U"? - )  +0 


3 


Now,  noting  that 


o„  U/2..  o 


qu(2 »K/fcr  ° 
K,  x  sin  e 

xo  0 


and  truncating  the  resulting  expression,  one  obtains: 


(C5) 


max 


xoKx  sine  0 

- 2_  +  |_  s1n(xo  +  l)e 

✓2ir  r  ° 


Along  e  =  */2  (C5)  reduces  to 


x  K. 

o  X 


max 


/2it  r 


o  (Kx  ^Ap 

T-  +  f-s1n(xo  +  l)f  =xo— V 

/2ir  r  0 


(C6) 


and  for  x  =  1/2,  the  LEFM  case,  we  recover  Eq.  B2. 

Since  one  expects  LEFM  to  prevail  away  from  the  free  surface,  one 

can,  by  assuming  x  =  1/2,  obtain  o°  from  the  slope  of  a  plot  of: 

o 

_(Kl>l%  VS  (I) 1/2  (see  Eq.  B4) 
o(ira) 

This,  in  effect  determines  the  value  of  the  non-singular  term  on 
R.S.  of  Eq.  (C6) . 

If  we  define  this  term  as  tq,  l.e. 

tq  =  f  s1n(xo  +  1)  f  (C7) 


then  Eq.  (C7)  suggests  that  tq  may  vary  If  xq  varies.  Thus,  as  one 
approaches  the  free  surface,  one  should  experience  a  variation  in  tq. 
LEFM  determinations  of  o°  show  that  it  begins  to  change  significantly  a 
short  distance  from  the  free  surface.  If  one  then  imposes  a  linear 

decrease  in  tq  from  this  point  to  zero  for  the  surface  slice,  one  finds 

* 

reasonable  correlation  between  the  photoelastic  and  moire  data  within 
the  transition  zone.  Using  this  approach,  and  Eq.  (C6)  in  the  form 


we  may  extract  values  of  \g  In  the  transition  region  from  a  plot  of 

109<W  -  T0>  vs  l09  r' 
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APPENDIX  D  CHOICE  OF  EXPERIMENTAL  METHODS 

Whole-field  optical  methods  potentially  applicable  for  these 
deformation  measurements  Include  the  moire, grid4  and  speckle® 
methods  for  In-plane  displacements;  and  shadow  moire, projection 
moire6  and  speckle5  methods  for  out-of-plane  displacements.  Durelli 
offers  additional  perspective  In  Ref.  7. 

Speckle  methods  have  been  applied  to  various  fracture  problems, 
notably  by  Chiang5  and  coworkers,  but  the  cases  treated  Involved  smaller 
displacements  than  encountered  here.  The  propellant  crack  problem 
Involves  displacements  that  are  too  large  to  maintain  speckle  correla¬ 
tion  between  the  before  and  after  deformation  Images,  so  speckle  tech¬ 
niques  are  Inappropriate. 

In-plane  moire  techniques  have  advantages  for  the  far-fleld  regions 
of  the  propellant-crack  specimen.  In  the  near-field,  however,  the 
rather  large  rotations  and  strains  would  cause  ambiguous  results.  The 
fact  that  the  large  rigid  Inclusions  create  local  displacement 
irregularities  would  further  confuse  the  moire  fringe  patterns. 

The  grid  method  appears  to  offer  the  best  features  for  In-plane 
measurements.  Although  more  effort  is  required  for  data  reduction,  the 


>v  > 


laboratory  work  for  acquisition  of  the  data  is  comparatively  easy.  The 


grid  method  can  be  applied  in  regions  of  very  large  extensions  and 


rotations  fully  as  well  as  regions  of  modest  deformation.  The  actual 


surface  displacements  are  sensed  Independently  of  Irregularities  caused 


by  the  rigid  inclusions.  (This  quality  is  distinct  from  moire,  which 


requires  smoothly  varying  deformations;  for  a  field  of  Irregular 


deformations,  fringe  contrast  is  diminished  and  ultimately  lost  as  the 


severity  of  irregularities  increases.)  The  grid  method  was  selected  for 


this  work.  Its  resolution  depends  upon  the  data  extraction  techniques 


employed  and  with  the  relatively  simple  techniques  used  here,  the 


resolution  was  estimated  (in  the  section  Results:  Global  Deformations) 


as  2.5  ym  or  100  yin.  Moire  with  about  2000  lines  per  inch  (80  lines 


per  mm)  would  be  required  for  equal  resolution,  and  use  of  such  fine 


gratings  would  be  exceedingly  difficult  for  these  tests. 


For  out-of-plane  measurements,  the  correlation  problem  again  pre¬ 


cluded  the  speckle  method.  The  shadow  moire  method  offers  all  the 


needed  features,  being  easy  to  employ  and  fitting  the  needed  sensitivity 


range.  It  is  Incompatible  with  the  grid  method,  however,  for  simul 


taneous  U,  V  and  W  data  collection  from  the  same  specimen  surface.  How¬ 


ever,  it  was  rationalized  in  the  section  Utility  of  W  Measurements  that 


W  need  not  be  measured  on  the  same  surface  as  U  and  V.  Consequently, 


shadow  moire  data  were  taken  for  many  tests  (Table  2)  on  the  back  sur¬ 


face  of  the  specimen. 


The  projection  moire  method  of  W  determination  is  compatible  with 


the  grid  method,  simultaneously  on  the  same  surface.  An  innovative 


method  was  developed  for  simultaneous  data  collection  in  a  single  photo¬ 


graph.  The  method,  described  in  Appendix  G,  is  somewhat  similar  to  that 


of  Chlang  and  Williams®.  However,  they  used  moire  for  U  and  V  Instead 
of  the  grid  method;  and  they  used  single-beam  projection  Instead  of 
double  beam,  with  consequent  reduction  of  sensitivity. 
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APPENDIX  E  THE  GRID  METHOD 
Printing  the  Grid 

For  printing  a  grid  on  the  surface  of  the  specimen,  the  basic  tool 

used  In  this  work  Is  a  commercially  available  metal  screen.  This  is  a 

thin  nickel  sheet  with  regularly  spaced  square  openings.  These  screens 

are  cataloged  under  the  name  Electroformed  Mesh  and  sold  by 

Buckbee-Meers 
245  East  6th  Street 
St.  Paul,  MN  55101 


The  screens  selected  for  this  work  were  grids  of  125  lines  per  inch  and 


500  lines  per  Inch  (5  and  20  lines  per  mm).  The  standard  nickel 
thickness  Is  0.0001  In.  (2.5  ym),  but  thicker  screens  of  0.0003  in. 

(7.5  ym)  were  obtained  by  special  order  for  this  work.  The  increased 
thickness  makes  handling  easier.  The  square  openings  were  wider  than 
the  vertical  and  horizontal  metal  bars  that  comprised  the  screen;  the 
screen  has  the  appearance  of  the  printed  grid  In  Fig.  13. 

Several  alternative  printing  processes  were  used.  In  all  of  them, 
the  screen  was  fixed  to  the  specimen  surface  to  act  as  a  mask  or 
stencil.  Reflective  material  was  deposited  through  the  openings  in  the 
screen.  When  the  screen  was  removed.  Its  Image  appeared  on  the  specimen 
as  bright  square  dots  on  the  dark  (nonreflective)  specimen  material. 

The  screen  could  not  be  reused. 

It  was  necessary  to  fix  the  screen  in  Intimate  contact  with  the 
specimen.  In  one  alternative,  the  screen  was  laid  on  the  tacky  surface 
of  the  specimen  and  pressed  down  lightly  with  a  roller.  Then  an  ultra 
thin  film  of  aluminum  was  deposited  on  the  surface  by  high  vacuum 
evaporation.  The  screen  was  peeled  off  to  complete  the  operation. 

An  equivalent  method  was  used  on  the  non-tacky  side.  In  this  case 
a  very  thin  film  of  liquid  epoxy  monomer  was  applied  to  the  specimen  to 
provide  a  sticky  surface.  Then  the  process  proceeded  as  before. 

While  this  printed  grid  gave  good  results,  some  Improvement  was 
desired.  Some  dark  blemishes  appeared  In  photographs  In  zones  where  the 
local  slope  of  the  specimen  was  great,  l.e.,  where  a  salt  crystal 
rotated  out-of-plane  as  the  specimen  deformed.  The  dark  blemish 
appeared  because  that  part  of  the  surface  did  not  reflect  light  Into  the 
camera  lens.  A  higher  degree  of  diffuse  reflection  was  desired  so  light 
would  reflect  In  a  wide  angle  of  directions. 


Titanium  dioxide  powder  has  strong  diffuse  reflection  and  it  was 
used  successfully.  In  this  case,  a  thin  layer  of  vacuum  grease 
(silicone  grease)  was  spread  on__the  specimen  surface  and  the  screen  was 
pressed  onto  the  grease.  T10  powder  was  sprinkled  onto  the  surface, 
pressed  Into  the  openings  using  a  ball  of  absorbent  cotton  and  excess 
powder  was  blown  off.  When  the  screen  was  peeled  off,  the  powder  was 
left  stuck  to  the  grease  In  neat  square  dots.  The  results  were 
excellent  for  125  line  per  Inch  (5  lines  per  mm)  grids,  but  the  method 
was  exceedingly  difficult  for  500  lines  per  inch  (20  lines  per  mm). 

Test  Apparatus  and  Procedure 

The  experimental  arrangement  Is  Illustrated  schematically  In  Fig. 
El.  The  specimen  with  a  grid  printed  on  it  was  Installed  In  a  tensile 
testing  machine.  It  was  stretched  at  a  constant,  predetermined  cross¬ 
head  speed;  the  machine  was  equipped  with  a  strip-chart  recorder  that 
graphed  the  load  vs.  crosshead  displacement  curve. 

The  grid  was  Illuminated  by  diffuse  white  light.  A  35  mm  camera 
was  critically  focused  on  the  grid.  The  camera  lens,  with  100  mm  focal 
length,  was  displaced  from  the  camera  body  using  extension  tubes  to 
obtain  approximately  1:1  magnification.  The  lens  aperture  was  stopped 
down  to  f/8,  but  the  effective  aperture  was  smaller  because  of  the 
extension  tubes.  Exposure  time  was  1/30  second,  using  Kodak  Technical 
Pan  film. 

The  camera  was  equipped  with  a  motor  drive  and  successive  photo¬ 
graphs  of  the  grid  were  taken  while  the  testing  machine  stretched  the 
specimen.  The  photographs  were  made  at  predetermined  time  Intervals, 


using  an  electronic  Interval  timer  to  trigger  each  exposure  and  advance 


Collimated 

Illumination 


the  film.  The  Interval  timer  also  sent  a  voltage  pulse  to  the  strip- 
chart  recorder;  this  produced  a  mark  on  the  load-deflection  curve  at  the 
instant  each  exposure  was  made.  The  raw  data  of  the  grid  method  Is 
Illustrated  by  Figs.  10  and  13. 

The  shadow  moire  tests  were  carried  out  simultaneously,  as 
described  In  Appendix  F. 

Data  Extraction 

Photographic  prints  of  the  deformed  grid  Images  were  produced  with 
approximately  10  times  magnification.  The  print  was  fixed  to  a 
digitizing  table  used  to  measure  the  x,y  coordinates  of  the  grid 
intersection  points.  The  points  chosen  are  indicated  In  Fig.  16a.  In 
this  operation,  a  magnifier  with  a  fine  cross-hair  was  positioned 
manually  over  the  grid  Intersection  point  and  a  button  was  pressed  to 
transmit  the  x,y  coordinates  of  the  cross-hair  to  a  computer  memory 
file.  The  digitizing  table  has  a  resolution  of  0.001  In.  (25  ym),  which 
translates  to  0.0001  In.  (2.5  ym)  resolution  on  the  specimen  scale. 

While  this  procedure  Is  tedious,  favorable  comments  can  be  made. 
After  Initial  experience  Is  gained,  the  time  required  to  digitize  one 
grid  pattern  was  about  one  hour.  The  overall  effort  was  small  compared 
to  that  needed  to  develop  an  automated  data  extraction  system.  In 
addition,  an  automated  system  would  have  lower  potential  for  accuracy 
and  reliability,  given  the  variation  of  line  shapes  near  the 
Intersection  points  and  the  blemishes  or  noise  In  the  photographs. 

Data  Reduction 

After  digitization,  the  data  for  each  grid  Intersection  point  was 
stored  In  computer  memory  files,  viz.,  the  column  and  row  numbers  that 


identify  the  point  and  the  x  and  y  coordinates  of  the  point.  Computer 
programs  were  developed  to  perform  several  operations  on  the  data. 
These  include 


•  Subtraction  of  coordinates  from  two  grid  systems,  to  obtain  the  U,V 
displacements  of  points  relative  to  the  zero  deformation  condition  (or 
other  previous  condition). 

•  Four-point  averaging  of  displacements. 

•  Deletion  of  lines  or  columns  of  data  to  Increase  gage  lengths  and 
affect  smoothing. 

•  Calculation  of  strains  e  ,  e  and  y  using  U,V  data  from  4  points. 

jr  x  xy 

•  Contouring,  to  connect  points  of  equal  displacement  or  strain  to  form 
contour  plots. 

The  results  are  Illustrated  in  Fig.  16  as  obtained  from  a 
printer. 

APPENDIX  F  SHADOW  MOIRE  METHOD 

The  technique  is  illustrated  on  the  right  side  of  Fig.  El.  A  rigid 
reference  grating  with  opaque  bars  and  transparent  spaces  is  fixed  next 
to  the  specimen.  The  oblique  Illumination  beam  casts  shadows  of  the 
opaque  bars  onto  the  specimen  surface,  thus  forming  a  grating  of  dark 
and  bright  lines  on  the  specimen.  Light  reflected  diffusely  from  the 
shadow  grating  passes  through  the  reference  grating  to  the  observer 
(camera).  Interaction  of  the  shadow  grating  and  reference  grating 
creates  moire  fringes,  fringes  that  depict  the  contour  map  of  the  gap 
between  the  warped  specimen  surface  and  plane  reference  grating  surface. 
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This  map  represents  the  relative  out-of-plane  displacement  of  each  point 
on  the  specimen. 


In  these  experiments,  T10  powder  was  rubbed  onto  the  specimen 
surface  to  Increase  Its  diffuse  reflectivity.  The  reference  grating  was 
held  against  the  specimen  by  a  thin  film  of  vacuum  grease.  The  grease 
was  applied  near  the  ends  of  the  reference  grating  and  squeezed  into  a 
thin  film  to  maintain  a  minimum  separation  between  the  grating  and 
specimen;  the  grease  did  not  extend  into  the  measurement  zone.  In 
addition  the  grating  was  restrained  from  vertical  translation  by  a  frame 
attached  to  the  lower  (fixed)  grip. 

The  relationship  defining  out-of -plane  displacements  W  is 


where  P  is  the  pitch  of  the  reference  grating,  f  is  its  frequency,  e  is 
the  angle  of  illumination  (Fig.  El)  and  N  is  the  fringe  order  in  the 
moire  pattern.  In  practice  W  must  be  interpreted  as  the  relative 
displacement  between  points  on  the  specimen,  since  there  is  no 
accounting  for  the  arbitrary  initial  gap. 

In  this  work,  f  was  20  lines  per  mm  and  0  was  45°.  The  sensitivity 
W/N  was  0.05  mm  (0.002  in.)  per  fringe  order.  An  example  of  the  shadow 
moire  fringe  pattern  is  given  in  Fig.  11. 


APPENDIX  G  U,V,W  METHOD 

This  method  captures  the  information  needed  for  whole-field  deter¬ 
mination  of  U,V  and  U  in  a  single  photograph.  As  before,  the  specimen 
is  prepared  with  a  125  lines  per  inch  (5  lines  per  mm)  grid,  using  the 
T10  powder  method.  Double-beam  projection  moire  is  used  to  create  a  low 


visibility  moire  pattern  of  W  on  the  specimen.  A  photograph  records 
this  pattern  on  a  background  of  the  deformed  grid,  as  shown  in  Fig.  61. 
The  U  and  V  displacement  information  can  be  extracted  exactly  as 
described  In  Appendix  E.  The  W  Information  Is  extracted  by  optical  fil¬ 
tering,  using  the  negative  from  which  Fig.  G1  was  printed. 

The  experimental  arrangement  is  illustrated  schematically  in  Fig. 

G2(a).  The  specimen  Is  Illuminated  by  symmetrical  beams  A  and  B.  These 

beams  do  not  have  constant  Intensity  across  their  widths,  but  Instead 

♦ 

they  exhibit  closely  spaced  dark  and  bright  stripes.  These  stripes  are 
projected  onto  the  diffusely  reflecting  specimen  and  they  are  photo¬ 
graphed  by  the  camera.  The  stripes  are  called  projected  gratings.  The 
two  projected  gratings  from  A  and  B  form  a  moire  pattern  of  W.  It  is  a 
moire  pattern  of  intensity  addition  (Gl),  however,  which  inherently  has 
poor  visibility  until  the  optical  filtering  process  is  used. 

The  method  used  to  form  the  projected  gratings  is  illustrated  in 
the  figure.  Here,  light  from  a  laser  is  expanded  into  a  conical  beam. 

It  passes  through  a  subsystem  comprised  of  the  beam-splitter  cube  BS  and 
mirrors  Ml  and  M2.  When  the  mirrors  are  parallel  to  their  adjacent  cube 
faces,  light  propagates  toward  the  parabolic  mirror  as  a  single  diverg¬ 
ing  beam.  However,  when  Ml  is  rotated,  two  separate  beams  emerge,  one 
shown  by  solid  lines  and  the  other  by  dashed  lines.  The  two  beams  are 
collimated  by  the  parabolic  mirror,  but  they  retain  their  angular  separ¬ 
ation  a.  Since  the  two  beams  are  mutually  coherent,  they  combine  in 
constructive  and  destructive  interference.  The  interference  fringes  are 
perpendicular  to  the  plane  containing  angle  a.  Their  frequency  (fringes 
per  unit  length)  is  proportional  to  a  when  a  is  small.  The  direction 
and  frequency  of  fringes  are  easily  adjusted  by  rotation  of  mirror  Ml. 


The  grating  frequency  that  appears  on  the  plane  of  the  specimen  depends 
also  upon  angle  e.  In  this  work  the  frequency  of  the  grating  projected 
on  the  specimen  was  200  lines  per  Inch  (8  lines  per  mm)  and  Its  orienta¬ 
tion  was  15°  from  the  y  axis,  as  Illustrated  In  Fig.  G2(b).  This  fre¬ 
quency  and  angle  were  obtained  by  adjusting  mirror  HI. 

Mirrors  M4  and  M5  are  used  to  project  the  second  grating  on  the 
specimen.  When  the  angle  of  B  Is  symmetrical  to  A,  the  frequencies  of 
the  two  projected  gratings  are  equal  and  a  null  field  Is  achieved.  The 
adjustment  Is  difficult,  however,  since  gratings  of  200  lines  per  Inch 
(8  lines  per  mm)  cannot  be  seen  with  the  naked  eye,  and  since  the  moire 
fringes  are  virtually  Invisible.  An  Indirect  technique  was  used  In  this 
work,  but  a  practical  means  of  adjusting  beam  B  to  achieve  an  Initial 
null  field  Is  to  view  the  moire  fringes  using  a  television  camera  and  a 
monitor  screen  adjusted  for  high  contrast.  The  specimen  should  be 
viewed  at  sufficiently  high  magnification  to  resolve  the  Individual 
lines  of  the  projected  gratings  on  the  monitor  screen.  Mirrors  M4  and 
M5  should  be  adjusted  In  this  operation,  not  Ml,  to  achieve  a  null  field 
prior  to  deforming  the  specimen. 

An  auxiliary  white  light  source  Is  used  to  take  the  photographic 
exposure.  This  enhances  the  visibility  of  the  grid  In  the  photograph. 
The  Intensity  of  the  auxiliary  light  was  adjusted  during  trial  exposures 
to  achieve  a  good  balance  between  the  grid  visibility  (for  U  and  V)  and 
the  projected  moire  Information. 

The  optical  filtering  arrangement  used  In  this  work  Is  Illustrated 
In  Fig.  G2(c).  The  negative  to  Fig.  G1  was  placed  In  the  object  plane; 
Its  enlarged  Image,  focused  by  the  parabolic  mirror,  appeared  In  the 
Image  plane.  A  camera  back  was  located  In  this  plane  to  record  the 
filtered  Image  on  photographic  film. 


The  negative  placed  In  the  object  plane  contained  an  Image  of  the 


projected  gratings.  This  Image  Is  Itself  a  grating,  capable  of 

diffracting  the  transmitted  light.  Light  from  each  of  Its  diffraction 

orders  appeared  In  the  filter  plane  as  bright  spots,  uniformly  spaced 

along  a  105°  line  (perpendicular  to  the  15°  projected  grating  lines). 

Light  from  the  +1  diffraction  order  was  allowed  to  pass  through  an 

aperture  (or  hole)  In  an  opaque  screen  placed  in  the  filter  plane.  All 

other  light  was  obstructed.  This  light  carried  the  projected  moire 

« 

pattern  (the  W  field),  with  enhanced  contrast,  to  the  film  plane.  The 
result  Is  Illustrated  by  Fig.  12.  This  Is  a  contour  map  of  the  out-of- 
plane  displacement  of  the  specimen  surface,  showing  Its  global 
deformation  plus  Irregularities  caused  by  local  elements  of  the 
heterogeneous  material. 

Light  was  also  diffracted  by  the  grid  lines  present  In  the  negative 
(object  plane).  These  diffractions  appeared  as  spots  In  the  filter 
plane,  but  the  spots  were  sufficiently  distant  from  the  aperture  to  be 
obstructed,  or  filtered,  out. 

The  basic  sensitivity  Is  twice  that  of  the  shadow  moire  method. 
Displacements  W  are  related  to  fringe  orders  N  by 
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where  f  Is  the  frequency  of  projected  grating  lines  on  the  specimen 
surface.  In  this  work,  f  was  200  lines  per  Inch  (8  lines  per  mm) 
and  e  was  65°,  to  yield  a  sensitivity  W/N  of  0.0012  Inches  per  fringe 
order  (0.03  mm  per  fringe  order). 
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APPENDIX  H  PHENOMENOLOGICAL  ANALYSIS  OF  CRACK  OPENING  AND  GROWTH 
Introduction 

Many  qualities  of  crack  opening  and  growth  cannot  be  observed 
directly,  but  can  be  Inferred  from  the  experimental  evidence.  A  process 
of  observations  and  self-consistent  deductions  was  undertaken  and 
reported  In  this  section.  Because  of  the  random  nature  of  the  material 
—  random  orientations  and  spaclngs  of  the  Inclusions  —  one  can  only 
generalize  on  Its  behavior.  The  Interpretation  Is  Intended  to  represent 
typical  behavior,  while  fully  realizing  that  diverse  variations  are 
present. 

Specimen  1  exhibits  the  characteristics  to  be  discussed.  Its 
deformation,  depicted  by  Fig.  13,  will  be  used  as  an  example. 


Crack  Blunting 

Figure  13a  shows  the  progress  of  crack  opening  and  blunting  as  the 
specimen  Is  stretched.  Image  (g)  shows  the  last  phase  of  blunting  on 
the  specimen  surface  prior  to  void  formation  and  crack  extension.  The 


material  adjacent  to  the  blunt  end  of  the  crack  suffers  very  large 
elongations.  The  zone  of  very  large  elongations  extends  ahead  of  the 
crack,  essentially  as  an  equilateral  triangle  with  the  blunt  crack  tip 
as  Its  base. 


The  deformation  of  the  surrounding  material  Is  qualitatively 
evident  from  observation  of  the  grid  and  quantitatively  from  Fig.  16. 
Above  and  below  the  crack  lips,  deformations  are  very  small.  Finite 


strains  appear  In  the  region  surrounding  the  triangular  zone,  where 
extensions  (e  )  and  shears  (y„v)  are  of  the  order  of  30*.  These  are 
average  strains  reflecting  the  Integrated  deformation  of  the  compliant 
rubber  and  stiff  (virtually  rigid)  Inclusions. 

An  explanation  or  model  consistent  with  the  experimental 
observations  Is  proposed  below. 

Localization  of  the  Pre-crack  Zone 

Figure  HI  Is  an  Idealized  schematic  view  of  the  deformation 
process.  The  original  crack  with  zero  opening  Is  shown  In  (a). 

Crystals  of  salt  are  Illustrated  as  squares  along  the  crack  boundary  and 
along  a  path  ahead  of  the  crack.  This  representation  Is  Idealized 
Inasmuch  as  the  crystals  are  given  an  orderly  alignment  and  spacing;  In 
addition,  the  secondary  system  of  rigid  Inclusions  are  not  Illustrated. 
The  salt  crystals  are  surrounded  by  the  rubber  matrix,  which  occupies 
all  the  space  not  occupied  by  the  crystals.  The  thickness  of  the  rubber 
Is  exaggerated,  too.  The  largest  crystals  In  the  Inert  propellant  are 
cubes  with  faces  more  than  0.4  mm  wide,  more  than  twice  the  pitch  of  the 
grid  In  Fig.  13,  whereas  the  thickness  of  the  rubber  Is  very  much 
smaller. 

As  the  specimen  Is  stretched  (perpendicular  to  the  crack  plane), 
the  rubber  Immediately  ahead  of  the  crack  suffers  the  largest  normal 
strain  ey.  The  basic  thesis  here  Is  Illustrated  In  Fig.  HI:  the  thin 
cross-hatched  zone  In  (a),  consisting  of  a  continuous  path  between 
crystals.  Is  extended  to  form  the  triangular  zone  shown  In  (b).  This 
will  be  called  the  pre-crack  zone,  since  the  crack  will  subsequently 
grow  Into  this  zone.  The  material  In  this  triangular  zone  Is  capable  of 
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sufficiently  large  extension  that  parent  rubber  outside  the  zone  Is  not 
required  to  share  the  burden  of  very  large  extensions. 

Figure  H2  Illustrates  another  Idealized  geometry.  In  which  the 
crystals  are  oriented  45°  to  the  crack  plane.  With  this  distribution  of 
crystals,  the  crack  Itself  must  have  jagged  lips,  essentially  as  shown. 
The  matrix  material  In  the  jagged  path  shown  by  crosshatching  spreads 
out  to  form  the  pre-crack  zone. 

The  orientations  and  spaclngs  of  Inclusions  In  the  real  specimen 
material  Is  far  more  random.  Nevertheless,  the  Inclusions  are 
surrounded  by  matrix  rubber,  and  continuous  paths  of  rubber  analogous  to 
those  Illustrated  here,  exist  for  any  configuration.  The  thesis  remains 
the  same:  the  trlagular  pre-crack  zone  Is  formed  from  rubber  that 
occupied  a  single  narrow  path  between  Inclusions. 

Rubber-rich  Surface  Layer 

The  specimen  material  used  here  was  cast  as  plates  of  uniform 
thickness  and  the  specimens  retained  the  original  as-cast  surfaces. 

These  surfaces  were  rich  In  rubber,  with  a  much  smaller  concentration  of 
salt  crystals  and  other  particles  than  the  Interior.  The  thickness  of 
the  surface  layer  appears  to  be  about  0.1  mm  on  the  tacky  surface  of 
Specimen  1  (see  Table  2).  The  rubber-rich  layer  was  thinner  on  the  non- 
tacky  surface  of  the  specimen  material,  approximately  1/4  to  1/2  as 
thick. 

Auxiliary  Observations 

When  the  specimen  material  Is  cut  or  fractured,  salt  crystals 
appear  Imbedded  In  the  fracture  surface.  They  can  be  pried  out  easily 


with  a  sharp  pick  and  the  crystals  usually  emerge  with  clean  surfaces. 

In  only  a  few  cases,  rubber  remains  attached  to  the  crystal  after  it  Is 
pried  out.  The  Implication  Is  that  the  bond  strength  between  the  rubber 
matrix  and  salt  crystals  Is  low  —  the  tensile  strength  and  shear 
strength  of  this  bond  Is  smaller  than  the  cohesive  strength  of  the 
rubber. 

Observations  were  made  with  a  bar  of  the  Inert  propellant  cut  from 
the  sheet  and  loaded  as  a  beam  In  bending.  The  maximum  tensile  stress 
was  developed  at  the  cut  surface;  the  bending  forces  were  parallel  to 
the  rubber-rich  surfaces  of  the  material.  The  beam  was  bent  and  viewed 
with  magnification  on  the  tensile  side.  At  a  sufficient  tensile  strain, 
separation  occurred  between  the  rubber  and  one  side  of  some  crystals. 
This  formed  voids  which.  In  the  case  of  the  larger  crystals,  were  easily 
seen.  The  voids  became  wider  In  the  direction  of  the  tensile  stress  as 
bending  proceeded.  Subsequently,  the  rubber  walls  between  neighboring 
voids  ruptured  and  a  cradk  was  formed  on  the  specimen  surface.  The 
crack  penetrated  Into  the  beam  and  at  some  depth,  narrow  ligaments  of 
rubber  could  be  seen  spanning  the  width  between  the  lips  of  the  crack. 
The  crack  grew  deeper  as  the  beam  was  bent,  and  this  growth  proceeded  by 
rupture  of  these  ligaments  and  creation  of  new  ligaments  (between  voids) 
at  greater  depths. 

The  rubber-rich  surface  layers  also  ruptured  as  the  crack  width  and 
depth  Increased.  However,  this  rupture  lagged  the  ruptures  of 
ligaments;  the  depth  of  the  crack  at  the  rubber-rich  surface  was 
typically  less  than  the  depth  at  which  Interior  ligaments  appeared.  The 
magnitude  of  the  lag  was  greater  on  the  tacky  thicker  rubber-rich 
surface  and  smaller  on  the  opposite  side  —  the  nontacky  thinner  rubber- 
rich  surface. 


Failure  Process 

The  foregoing  observations  used  flexural  specimens  so  that  the 
failure  process  could  be  viewed  from  a  direction  perpendicular  to  the 
crack  front.  While  the  specimen  geometry  was  different,  the  basic 
features  of  failure  are  believed  to  be  the  same  as  those  In  the  biaxial 
tension  tests. 

Figure  H3  Illustrates  the  hypothesized  fracture  process  In  the 
biaxial  specimen  of  Fig.  13,  for  a  highly  Idealized  crystal  arrangement. 
In  (b),  the  pre-crack  zone  Is  characterized  by  large  voids  and  narrow 
ligaments  separating  the  crystals.  The  matrix  rubber  Is  separated  from 
opposing  horizontal  faces  of  the  crystals  to  form  elongated  box-shaped 
voids.  The  rubber  Is  also  unbonded  from  the  vertical  sides  of  opposing 
crystals  and  this  material  Is  stretched  to  form  the  ligaments  and  the 
extended  surface  layer.  The  voids  are  closed  cells,  l.e.,  the  ligaments 
are  hollow  tubes  of  matrix  rubber.  In  (c)  the  crack  opening  Is  great 
enough  to  cause  rupture  of  ligaments  but  the  stretched  surface  layer  Is 
Intact.  With  greater  crack  opening  (d)  the  surface  layer  Is  ruptured 
and  the  rubber  relaxes  back  to  approximately  Its  original  shape. 

Figure  H4  Illustrates  an  Internal  section  of  the  specimen.  The 
separate  conditions  of  Fig.  H3  occur  at  different  regions  of  the 
specimen.  Identified  here  by  corresponding  letters  In  parentheses.  The 
ligaments  appear  as  very  narrow  zones,  but  a  three-dimensional  view 
would  show  them  as  thin  walls  of  elongated  tubes. 

In  Fig.  H4,  notice  the  two  tiny  crystals  located  between  two 
vertical  elements  of  rubber.  For  such  a  case,  stretching  of  the 
surrounding  rubber  matrix  would  carry  the  lower  crystal  down  Into  the 
triangular  pre-crack  zone.  The  presence  of  crystals  In  the  ligament 
zone  was  not  detected,  however.  In  microscopic  Inspections. 
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An  Important  question  arises:  Why  Is  the  crack  end  so  blunt?  If 
the  shear  strength  between  the  rubber  and  crystals  Is  low,  why  doesn't 
the  rubber-rich  surface  layer  act  Independently  as  a  thin  sheet  of 
rubber,  which  would  form  an  elliptical-shaped  crack  opening?  The  answer 
lies  in  Fig.  H3(c).  While  the  surface  layer  debonds  from  the  vertical 
sides  of  the  opposing  crystals.  It  remains  connected  to  the  horizontal 
elements  of  rubber  matrix.  This  keys  the  surface  layer  to  the  crystals 
on  both  sides  of  the  pre-crack  zone  and  it  stretches  between  these 
points  to  form  the  nearly  straight  crack  tip. 

Referring  to  Fig.  H3  (b  and  c),  take  note  of  the  triangular  por¬ 
tions  of  the  rubber-rich  surface  layer.  Strains  In  the  surface  layer 
can  change  rather  abruptly  near  the  crystals,  where  the  rubber  Is  keyed. 
Elsewhere  In  the  surface  layer,  the  strains  must  change  In  a  smooth 
transition  from  the  smaller  strain  In  the  cross-hatched  region  to  the 
larger  strain  in  the  vertically-hatched  region.  This  transition  occurs 
In  the  triangular  portions  of  the  rubber-rich  layer. 

Some  numerical  values  will  help  give  perspective  to  these  condi¬ 
tions.  Crystals  exceeding  400  ym  in  width  were  found  In  abundance  In 
the  material.  Assume  a  representative  size  of  all  crystals  Is  100  ym 
and  assume  a  uniform  array  of  100  ym  cubical  crystals,  each  crystal  sur¬ 
rounded  uniformly  by  the  rubber  matrix.  For  75%  volume  fraction  of 
crystals,  the  thickness  of  the  rubber  matrix  Is  calculated  as  10  ym. 

The  cross-sectional  areas  of  rubber  and  crystal  are  In  the  ratio  1:5,  or 
1/6  the  total  area  Is  rubber. 

The  thickness  of  the  rubber-rich  surface  layer  (on  the  tacky  side) 
Is  about  equal  to  the  representative  cube  width.  The  crack  blunting  In 
Fig.  13a  (Image  g)  Is  about  1  mm.  If  the  conditions  of  Fig.  H3(b) 


prevail,  the  extended  length  of  the  surface  layer  is  1  mm  plus  two 
crystal  heights,  or  1.2  mm;  Its  original  length  [Fig.  H3(a)l  Is  0.2  mm. 
The  simple  tensile  strain  In  the  surface  layer  at  the  blunt  end  Is  ey  ■ 

5  or  500£.  This  Is  reasonable. 

A  similar  model  of  voids  and  connecting  ligaments  Is  described  In 
Ref.  HI  by  Bills  and  Blatz.  They  remark,  however,  that  certain  data 
challenges  the  theory.  Whereas  the  tensile  strength  of  the  material 
would  be  expected  to  decrease  with  the  volume  fraction  of  Inclusions, 
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they  found  essentially  constant  strengths  of  HTPB  propellant.  Indepen¬ 
dent  of  the  volume  fraction.  This  factor  remains  unexplained.  Similar 
tensile  tests  were  not  conducted  for  the  present  Inert  propellant. 

Force  Transfer 

Different  mechanisms  of  force  transfer  should  exist  In  the  trian¬ 
gular  pre-crack  zone  (Including  regions  Immediately  surrounding  this 
zone)  and  the  remaining  regions.  Remote  from  the  high  strain  regions, 
one  could  reasonably  hypothesize  that  debonding  between  crystals  and 
matrix  rubber  does  not  occur.  The  condition  Is  modeled.  Ideally,  by 
Fig.  H3(a).  Under  a  vertical  tensile  force,  the  vertical  matrix 
elements  are  restrained  from  stretching  by  the  adjacent  rigid  crystals. 
The  horizontal  matrix  elements  become  thicker,  or  exhibit  a  local  strain 
ey.  This  accounts  for  some  small  average  strain  In  the  composite 
material. 

At  a  higher  strain  level,  bonds  fall  at  the  horizontal  crystal- 
matrix  Interfaces.  This  Is  because  of  strong  trlaxlal  tensile  stresses 
developed  adjacent  to  the  rigid  crystals,  as  In  the  classical  poker-chip 
test;  the  rubber  cannot  respond  with  a  change  of  Its  volume  and  failure 


At  still  higher  strain  levels,  debonding  occurs  by  shear  at  the 
vertical  Interfaces  between  crystals  and  matrix.  Subsequent  strain  Is 
accommodated  by  stretching  of  these  vertical  rubber  elements,  while  the 
horizontal  elements  remain  unloaded. 

A  progression  is  prescribed  In  which  load  Is  carried  exclusively  by 
horizontal  matrix  elements  at  first  and  later  exclusively  by  vertical 
elements.  Strains  vary  accordingly.  In  low  strain  regions,  the 
vertical  tensile  load  Is  carried  on  a  continuous  horizontal  path  of 
rubber  of  area  A.  In  high  strain  regions  the  load  Is  carried  by  the 
vertical  matrix  elements  only,  whose  representative  area  Is  estimated  In 
the  previous  section  as  A/6.  Tensile  stress  accumulates  with  tensile 
load  six  times  faster  In  the  pre-crack  zone  than  In  the  low  strain  zone. 
The  highly  nonlinear  load-displacement  curve  of  Fig.  10  reflects  this 
quality. 

What  condition  exists  Immediately  adjacent  to  the  triangular  pre¬ 
crack  zone,  where  the  average  strains  are  of  the  order  of  30%?  With  the 
representative  10:1  crystal  to  matrix  widths,  an  assumption  that  the 
strain  Is  concentrated  In  the  horizontal  rubber  elements  of  Fig.  H3(a) 

Is  unreasonable;  since  the  crystals  are  rigid,  this  would  lead  to 
strain  e  In  the  rubber  exceeding  300%.  Instead,  the  strain  must 
represent  30%  change  of  length  of  the  vertical  matrix  elements.  This 
requires  debonding  from  the  rigid  crystals  and  also  void  formation  along 
the  horizontal  crystal  faces.  The  condition  would  be  represented  by 
Fig.  H3(b)  where  the  voids  are  much  smaller,  viz.,  void  height  about  30% 
of  crystal  height.  In  this  zone,  like  the  pre-crack  zone,  force 
transfer  Is  through  the  vertical  matrix  elements. 


While  these  generalizations  are  developed  for  an  orderly  array  of 
crystals.  It  appears  reasonable  that  they  apply  In  principle  to  random 


arrays  and  to  the  specimen  material. 

Crack  Growth 

Figure  13b  shows  a  long  sequence  of  crack  growth  for  the  same 
specimen  as  Fig.  13a.  The  fixed  parameter  was  cross-head  speed,  or  rate 
of  separation  of  the  specimen  grips.  The  load  on  the  specimen  was 
decreasing  monotonlcally  as  shown  In  Fig.  10. 

In  Fig.  13b  [Images  (a)-(d)],  two  voids  appear  In  the  pre-crack 
zone,  but  they  subsequently  become  part  of  the  main  crack  by  disappear¬ 
ance  of  material  to  their  left.  In  Images  (e)-(h),  the  pre-crack  zone 
advances  along  an  oblique  path  and  voids  are  formed  which  subsequently 
coalesce  with  the  main  crack.  In  Images  (1)-(1),  a  new  pre-craze  zone 
Is  formed  at  a  small  distance  from  the  main  crack;  voids  form  and 
enlarge,  but  they  are  separated  from  the  main  crack  by  a  bridge.  The 
bridge  continues  to  stretch  until  Image  (m);  by  that  stage  the  bridge 
material  has  failed  and  disappeared,  and  a  new  pre-crack  zone  has  formed 
at  the  blunt  end  of  the  crack.  In  Images  (n)-(p),  we  see  another 
episode  of  crack  extension  and  void  formation.  The  dramatic  qualities 
revealed  by  the  sequence  of  photos  are  huge  shear  strains,  void  forma¬ 
tion  and  coalescence,  path  changes  and  bridging.  These  will  be  ration¬ 
alized  In  the  following  sections. 


Oblique  Paths  and  Shear  Strains 

The  pre-crack  and  subsequent  crack  can  extend  along  an  oblique 
path,  as  seen  In  Fig.  13b,  Images  (e)-(h).  Figure  H5  Illustrates  the 
mechanism  consistent  with  the  observations  and  the  stated  thesis. 
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Again*  the  pre-crack  zone  Is  created  by  a  vertical  translation  that 
separates  the  upper  and  lower  boundaries  of  the  vertically  translated 
crystals.  A  horizontal  and  a  vertical  grid  line  are  shown  crossing  the 
crack  path  In  (a)  and  the  pre-crack  zone  In  (b).  The  original  positions 
of  points  A  and  B  are  changed  by  a  relative  displacement  equal  to  the 
vertical  translation;  similarly,  C  and  D  are  displaced  vertically.  One 
would  expect  to  see  sharp  discontinuities  In  the  horizontal  grid  line, 
as  Illustrated  by  the  full  line  In  (b).  This  represents  the  deformation 
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presumed  to  exist  In  the  crystals. 

The  surface  deformation  Is  different,  however,  as  a  result  of  the 
rubber-rich  surface  layer.  Referring  again  to  Fig.  H3(b  and  c),  large 
deformations  on  the  outside  surface  do  not  commence  opposite  the 
translated  crystal  boundaries,  but  at  points  above  and  below  the 
boundaries.  This  means  that  large  vertical  displacements  commence  at 
some  point  A'  In  Fig.  H5b,  a  point  outside  the  pre-crack  boundary.  The 
vertical  displacement  builds  upbore  gradually  on  the  outside  surface; 
consequently,  the  originally  horizontal  grid  line  dips  down  gradually 
along  path  A'O,  and  It  recovers  Its  horizontal  slope  gradually  along 
path  OB'. 

This  dramatic  curvature  of  the  Initially  horizontal  lines  Is  seen 
In  Fig.  13b,  Image  (e)-(h).  Their  slopes  become  nearly  vertical,  nearly 
parallel  to  the  Initially  vertical  grid  lines.  This  signifies  severe 
shear  strains,  which  exist  In  the  rubber-rich  surface  layer. 

Failure  of  Ligaments  and  Surface  Layer 

The  crack  growth  Is  explained  as  a  progression  of  ligament 
extensions  and  failures.  On  the  outside  surface,  however.  It  Is  a 
progression  of  void  formation  and  coalescence  In  the  rubber-rich  layer. 


Failure  of  the  ligaments  could  be  expected  to  occur  at  a  limiting 
strain  (assuming  uniform  strain  rate).  Typically,  the  limiting  strain 
would  be  reached  first  at  the  blunt  end  of  the  triangular  pre-crack 
zone.  Exceptions  are  possible,  however,  because  the  original  lengths  of 
the  ligaments  are  not  all  the  same.  The  original  length  relates  to  the 
vertical  debond  length  between  crystals  and  rubber,  beyond  the  position 
where  the  rubber  Is  keyed  and  fixed  by  crystals;  see  Figs.  H3  and  H4. 
With  random  crystal  arrangements,  this  debond  length  must  be  random, 
too,  causing  the  strain  to  vary  from  ligament- to- ligament.  Conse¬ 
quently,  ligament  failures  will  not  advance  exclusively  from  the  blunt 
end  of  the  crack;  some  will  occur  further  ahead  of  the  blunt  end. 

Failure  of  the  rubber-rich  surface  layer  Is  analogous.  Figure  13b 
shows  that  voids  typically  appear  where  tensile  strains  and  shear 
strains  are  both  large.  Again,  the  magnitude  of  these  strains  Is 
related  to  the  local  conditions,  especially  the  positions  where  the  sur¬ 
face  layer  Is  keyed  to  crystals.  With  random  crystal  distribution,  the 
limiting  value  of  combined  tensile  and  shear  strain  can  occur  within  the 
triangular  pre-crack  zone  as  well  as  at  the  blunt  end  of  the  zone. 

What  Is  this  limiting  value?  Why  are  the  voids  elliptical  with 
major  axes  parallel  to  the  tensile  direction?  Probable  answers  were 
found  from  crude  tests  using  broad  rubber  bands.  Tensile  failure  sur¬ 
faces  were  perpendicular  to  the  load,  suggesting  that  maximum  principal 
stress  or  strain  Is  the  failure  criterion.  Figure  H6  Illustrates  that  a 
shear  strain,  when  combined  with  a  large  normal  strain,  significantly 
Increases  the  principal  strain,  but  changes  the  principal  strain  direc¬ 
tion  only  a  small  amount.  Consequently,  failure  leading  to  a  void  In 
the  rubber-rich  surface  layer  would  be  a  localized  rupture  or  crack  with 


a)  Grid  lines  for  the  case  of  c.  ■  3  and  c 


(b)  same  condition  except  Is  large.  The  bold  dashed  lines 


represent  maximum  extensions  and  maximum  principal  strain 


resents  a  crack  in  a  rubber  sheet  an 


represents  the  sheet  when  Its  height  is  extended  by  an 


strain  *4. The  crack  lips  deform  into  a  ta 


elliptical  boundar 


-  \  »■  ,  >  '  V  , 


•V-Ts’V'v' .*  •'  .  '•  •  *'  -vi- 


the  crack  plane  perpendicular  to  the  principal  strain.  This  Is  essen¬ 
tially  perpendicular  to  the  direction  of  tensile  loading. 

When  such  a  crack  Is  cut  In  a  broad  rubber  band  and  the  band  Is 
stretched,  the  crack  opens  Into  an  ellipse.  For  large  strains,  the 
ratio  of  lengths  of  the  vertical  (the  load  direction)  and  horizontal 
axes  of  the  ellipse  is  approximately  equal  to  the  principal  strain 
ratio.  This  is  Illustrated  In  Fig.  H6(d),  where  e  »  3.  Correspond- 
Ingly,  the  elliptical  voids  In  Fig.  13b  exhibit  nearly  vertical  prlncl- 
pal  axes. 

This  model  of  void  formation  and  coalescence  In  the  specimen  sur¬ 
face  Is  one  In  which  small  cracks  are  developed  In  the  rubber-rich  sur¬ 
face  layer  ahead  of  the  main  crack.  In  the  unloaded  condition  the  crack 
planes  would  be  perpendicular  to  the  maximum  principal  strain  direction 
—  essentially  the  y  direction,  the  loading  direction.  In  the  highly 
strained  condition,  the  lips  of  each  small  crack  extend  Into  an  ellip¬ 
tical  shape  with  major  axes  essentially  parallel  to  the  y  direction. 
These  cracks  grow  longer  (In  the  x  direction)  as  the  specimen  Is 
stretched,  causing  the  elliptical  void  to  enlarge.  When  adjacent  cracks 
join,  their  corresponding  voids  coalesce. 

Path  of  Crack  Growth 

Let  Fig.  H7  represent  (an  Idealized  schematic  diagram  of)  a  blunt 
crack  and  pre-crack  zone,  leading  to  the  candidate  path  of  crack  exten¬ 
sion.  How  can  we  rationalize  the  probable  path  of  the  crack  as  It 
extends  among  the  crystals? 

Returning  to  the  reasoning  of  Fig.  H3,  the  onset  of  debonding  along 
vertical  crystal /rubber  Interfaces  follows  debonding  at  horizontal 
Interfaces.  The  global  strain  level  at  which  horizontal  debonding 
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occurs  Is  Inversely  related  to  the  thickness  of  the  horizontal  elements 
of  matrix  rubber.  Accordingly,  the  onset  of  debonding  along  vertical 
Interfaces  Is  hastened  In  regions  where  the  horizontal  portions  of 
rubber  are  thinner  than  elsewhere. 

Where  debonding  along  vertical  surfaces  commences,  the  local 
stiffness  of  the  composite  Is  greatly  reduced  and  further  deformation 
concentrates  In  that  neighborhood.  By  this  argument,  the  path  of  the 
pre-crack  zone,  and  subsequently  the  path  of  the  crack  Itself,  follows 
the  route  along  which  the  vertical  height  of  matrix  material  Is 
smallest.  In  the  case  of  Fig.  H7,  the  upper  branch  beyond  point  A  would 
be  the  path  of  crack  extension. 

Crack  Stoppers,  Bridges 

Another  fascinating  behavior  Is  observed  In  the  deformation 
sequence  seen  In  Fig.  13b  (Images  (j)-(m)l,  deformation  that  proceeds  as 
the  specimen  Is  stretched.  In  Image  (j),  the  main  crack  has  achieved  a 
certain  length,  but  subsequent  views  show  no  Increase  of  Its  length. 
Instead,  a  new  pre-crack  zone  opened  slightly  below  the  tip  of  the  main 
crack.  Voids  grew  and  coalesced  In  that  zone  until  It  was  separated 
from  the  main  crack  by  a  bridge  of  highly  strained  material.  By  (m)  the 
bridge  failed  and  the  new  zone  joined  the  main  crack. 

An  explanation  of  crack  arrest  Is  required.  The  crack  extends 
Inexorably  along  a  single  path,  but  the  path  can  lead  to  a  crack  stopper 
—  like  a  dead-end  street.  This  Is  the  condition  where  a  local  crystal 
configuration  Inhibits  the  development  of  large  strains.  Such  a 
condition  Is  Illustrated  schematically  In  Fig.  H8.  Meanwhile,  a  point  B 
will  exist  In  the  neighborhood  of  A,  where  local  strains  are  large,  and 


a  new  crack  path  Is  established.  The  localization  of  deformation 
continues  along  this  new  path  until  It  leads,  by  chance,  to  a  crack 
stopper  —  to  a  local  zone  where  large  strains  cannot  occur. 
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